Ehrenfest time in the weak dynamical localization 
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The quantum kicked rotor (QKR) is known to exhibit dynamical localization in the space of its 
angular momentum. The present paper is devoted to the systematic first-principal (without a regu- 
larizer) diagrammatic calculations of the weak-localization corrections for QKR. Our particular em- 
phasis is on the Ehrenfest time regime - the phenomena characteristic for the classical-to-quantum 
crossover of classically chaotic systems. 
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I. INTRODUCTION 

In recent years it became abundantly clear that driven 
quantum systems exhibit behavior that is qualitatively 
different from their classical counterparts. Indeed, a 
driven classically chaotic system exhibits Joel's heating. 
That is, its average stored energy increases at a steady 
rate. In other words, such behavior may be character- 
ized as a diffusion in the system's phase space. The 
remarkable feature of driven quantum systems is finite- 
ness of their phase space motion (localization) . Such dy- 
namical localization phenomena was discussed in context 
of pumped quantum dotg^, Bose-Einstein condensates 
subject to pulses of optical standing wav o^i'^i'* , optical 
microcavity^, and other systems. 

The simplest model that became a paradigm for stud- 
ies of the quantum dynamical localization is quantum 
kicked rotor (QKR). It was numerical discovery of local- 
ization in QKR by Casati et al&^ in the late seventies 
that triggered the broad interest to the subject. Recent 
progress in trapping of cold atoms and optical manipu- 
lation with them lead to experimental realization of the 
QKR with the unprecedented degree ofcontrol^-l The 
kicked rotor is described by the time-dependent Hamil- 
tonian: 
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where angle 9 and angular momentum I are the pair of 
canonically conjugated variables. The amplitude of the 
kicks is described by the dimensionless parameter K , also 
known as the classical stochastically parameter. It is the 
only parameter of the corresponding classical problem. 
The quantum problem possesses another dimensionless 
parameter: the effective Planck constant k. The latter 
enters the problem through the canonical commutation 
relation: [6, 1] = ik. The two parameters K and k are 
straightforwardly related to the optical wavelength, am- 
plitude and atomic mass in cold atoms experiments^*^. 

Historically, the classical kicked rotor, or stan- 
dard mapping, first introduced by Chirikov, served as 
the prototype model for various transport processes in 



plasma o?!^? . It was established that the classical dy- 
namics of the kicked rotor exhibits complicated behav- 
ior. Most notably it demonstrates the transition from 
the regular to chaotic motionSiiS. In particular, for suffi- 
ciently large classical parameter K 5), the chaotic dif- 
fusion takes place in the space of angular momentumSiiS. 
The latter is associated with the unlimited diffusive ex- 
pansion of an initially narrow momenta distribution: 
S{P{t)) = 2Dcit . The classical diffusion constant, Dd, 
was a subject of numerous studies^' and is well under- 
stood by no'Viiiiii^. For large stochastically parameter 
if > 1, one finds D^i ~ + 0{K^/^), where the 

omitted corrections posses an oscillatory dependence on 
K. 

The pioneering numerical studies of Casati et al^^ re- 
vealed that the corresponding quantum system, S ^ 0, 
behaves in a dramatically different way. The initial diffu- 
sive expansion (that is heating) saturates after a certain 
time, t]^, and the momentum distribution width stabi- 
lizes at 6{P{t)) ~ = DciIl- It was soon suggested 
in Ref. that similarly to Anderson localization, quan- 
tum phase interference may lead to the suppression of 
classical diffusion for long enough time. This heuris- 
tic idea was complemented by mapping QKR onto a 
one-dimensional tight-binding Anderson model with the 
pseudo-random potentiaU^*i^. Such interpretation leads 
to the estimate of the localization length as ^ — Dd/k 
and thus ti, — Dd/k^- The similarity was further con- 
firmed by studies of a perturbation that breaks the "time- 
reversal symmetry" (TRS) of the QKlU^. Such pertur- 
bation suppresses the survival probability by a factor of 
two. That is closely analogous to the doubling of the 
Anderson localization length by the static magnetic field, 
destroying the interference between a trajectory and its 
time-reversal partne r^^i^° . 

If the momentum localization length, ^, is much longer 
than the "microscopic" scale of the angular momentum 
(that is k) then K ^ k and thus ^ 1. In this case 
there is a parametrically long crossover regime from the 
classical diffusion: S{P{t)) « 2Dcit for 1 < t <C to 
the strong localization: S{P{t)) — for t > t^. One 
may be able to develop a systematic perturbation the- 
ory in powers of (i/ti) <C 1, analogous to the weak- 
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localization loop expansion in the Anderson localization 
theorjiii. Such task was undertaken by Altland^^, who 
found for the one-loop correction: S{P{t)) = 2Dcit{l — 
0.75y/t/tL )• It was suggested furthermore that the uni- 
versal long-time behavior of the QKR is described by the 
diffusive supersymmetric nonlinear cr-modeliS similar to 
those employed in the localization theory^". 

Those calculations essentially map the QKR on a quan- 
tum particle in the field of a white-noise random poten- 
tial. While such analogy is reasonable at long time scales 
i ~ ^ 1, it fails to recognize details of the classical 
to quantum crossover at intermediate time scales. In- 
deed, a quantized classically chaotic system requires a 
certain time scale, called the Ehrenfest time, te, to de- 
velop quantum interference effects. This fact was real- 
ized independently in various context a^n-^i^ and now days 
is well documented in the literaturo'^'i^^i^^i^^i^^i^^i^'''. The 
physics behind this fact is as follows. To experience the 
quantum interference, two classical trajectories must con- 
verge to a region of the phase space of the size S166 « k. 
Convergence (and divergence) of trajectories in a classi- 
cally chaotic system is governed by the Lyapunov insta- 
bility exponent A as e.g. S9{t) ^ exp{— At}. It thus takes 
time tE ~ A~^ ln(l/^) before the interference effects can 
reveal themselves. 

As we show below, for iiT ^ 1 the quantitative def- 
inition of the Ehrenfest time for the QKR is given by 



tE 
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while the classical Lyapunov exponent is A = \n(K/2)^. 
Therefore there is the parametric regime 1 <C <C 
or more precisely: 



ln( - ) >lnJ-s:>0, 



(3) 



where exists a wide separation of the relevant time scales: 



K ifi < ti . 



(4) 



One may expect that such regime is amenable for an an- 
alytical treatment of the classical-to-quantum crossover. 
This problem was first tackled by Aleiner and Larkin^- 
in the context of random classical (long-range) potential 
scattering (e.g., random Lorentz gases). However, due 
to the complexity of the Lorentz gas classical dynamics, 
their treatment required a regularization. The latter is 
essentially a weak quantum scattering potential added to 
the Lorentz gas. 

The purpose of this paper is to develop a systematic 
first-principle analytic treatment of the QKR. In partic- 
ular, we are able to incorporate the semiclassical dynam- 
ics at the scale of t^; into the weak dynamical localiza- 
tion theory without introducing any regularization. The 
QKR allows thus to demonstrate explicitly an essential 



point: existence of the dynamical localization is an intrin- 
sic property of quantized classically chaotic systems - not 
an artefact of an extraneous regularization. (Remarkably, 
the Ehrenfest time does not dependend on the regular- 
izer strength2SiS2i2£.) This observation is fully consistent 
with the early studies of Ehrenfest timeSliSS, suggested 
the existence of such time scale. For the time interval 
t t^, our approach fully encompasses the Ehrenfest 
regime. The results were reported in the short letter—. 
The main result for the time-dependent spread of the 
wave-packet may be formulated as: 
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where d[t) is the standard Heaviside step-function. At 
intermediate times, tE ^ t < tL, Eq. lO crosses over to 
the standard weak-localization correctioni^. At shorter 
times, t ~ tE, there is a delay in developing localization 
given by AtE- A few comments are in due: (i) the actual 
delay is not absolutely sharp, as suggested by Eq. |(SJ). 
There are exponentially small deviation from the straight 
line 2Dcit even for t < AtE, which exact shape is calcu- 
lated below, (ii) Eq. jSJl describes quantum correction 
linear in k that appears to be delayed by AtE- As first 
noticed by ShepelanskjsSi, quadratic in k corrections show 
up even at earlier time. However, for at least three first 
kicks, they may be fully absorbed into a renormalization 
of the diffusion constanli^ {SDd ~ k'^, essentially due to 
the change in the scattering cross-section). It is thus an 
oversimplification to claim the absence of quantum effects 
at t < AtE- (in) Eq. ||SJ) constitutes the one-loop weak- 
localization correction. Below we report also the results 
of the two-loop calculation. It brings the next order cor- 
rection ~ 6{t — QtE){t — QtE)'^/tL that also "protects" 
the early time evolution from the localization effects. It 
is still an open problem to sum up the entire series to 
develop a theory of strong localization that accounts for 
the Ehrenfest time phenomena. 

Most of the existing experiments on atomic gases^i^SiSi 
do not fall down in the parametric regime © , but rather 
have ^ « 1. In this case « 1 and our result, Eq. ||SJ), 
can only be viewed as a qualitative one. We discuss below 
other possible realizations of the QKR that utilizes driven 
Josephson junctiona^^i^. Such systems may prove to be 
more suitable for exploring the parametric regime Q and 
thus for a quantitative comparison with the theory. 

The outline of the rest of this paper is as follows: 
Sec. m is devoted to a qualitative semiclassical picture 
of the weak-dynamical localizations and the Lyapunov 
regime. Sections IIIII IIVI and serve to quantify these 
ideas. In Sec. IIIII the diffusion in the phase space of 
the kicked rotor is obtained as a classical approximation 
to the full quantum propagation. Sec. IIVI is the cen- 
tral part of the present work. It formulates a general 
framework to deal with the weak dynamical localization 
at the semiclassical level. In particular, we calculate the 
frequency-dependent one-loop correction to the classical 
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diffusion coefficient and study its effect on the momen- 
tum dispersion. This formahsm is applied in Sec.[V]to a 
modified QKR with broken time-reversal symmetry. The 
frequency-dependent quantum corrections are calculated 
at the two-loop level. Experimental realizations of some 
driven quantum systems are discussed in Sec. IVII The 
effects of noise and dephasing are subject of Sec. IVIII 
We conclude in Sec lVIIll Some technical details are del- 
egated to Appendices. 



II. QUALITATIVE CONSIDERATIONS 

The physics of the weak-localization corrections is tra- 
ditionally discussed in the language of classical trajecto- 
ries. The classical motion of a particle in the random 
potential is characterized by a rapid randomization of 
momenta and diffusion spreading of the coordinate. It is 
thus customary to visualize a trajectory in the coordinate 
space as a random motion between static impurities. It 
is straightforward to develop a similar approach for the 
QKR. In the kicked rotor problem the roles of coordinate 
and (angular) momentum are interchanged. Indeed, for 
K ^ 1 the angular coordinate 9 is rapidly randomized 
(over the interval [— tt, tt]), while the angular momentum I 
acquires a (quasi) random change £ [—K,K]. The latter 
results in the diffusion in the space of angular momentum 
(see below). We shall thus visualize a "trajectory" as a 
sequence of values of the angular momentum that kicked 
rotor "visits" upon successive kicks, Fig. 

For a quantitative description of the classical motion 
it is convenient to monitor pairs of angle and angular 
momentum in discrete moments of time. This maps the 
classical dynamics onto a so-called standard map: 

In+l = In + Ksm9n, 

9n+l = 9n + ln+1 ■ (6) 

Notice that Z„ stands for the angular momentum immedi- 
ately after (n— l)-th kick, and 0„ for the angle before n-th 
kick. It is now obvious that the two successive points of 
the trajectory: l„ and In+i differ by K sin 9n- As a first 
approximation one may treat 6'„ as uniformly distributed 
over [0, 2tt] and thus ((Z„+i - /„)2) = K^isi-a^ 9) = K'^ /2. 
As a result, ((/„ - Zq)^) = 2Dcin with Dd = An 
account of the residual correlations between successive 
0„'s, leads to a renormalization of Dd with the next term 
scaling as K^/"^ etciiS*ii. For completeness of the presen- 
tation we shall derive the full result for the classical diffu- 
sion constant in Appendix^ For the current qualitative 
discussion it is enough to appreciate that a trajectory 
of the classical kicked rotor exhibits random hops in the 
space of angular momentum, leading to: 

5{l^{t)) ^2Ddt. (7) 

This result has a simple physical interpretation: the av- 
erage energy of the kicked rotor linearly increases with 



time, reflecting a constant rate Joel's heating. It is ex- 
actly this property of the classical kicked rotor that made 
it useful in the accelerator physics^iili. 




FIG. 1; The first quantum correction to the density-density 
correlator: (a) one-loop weak localization diagram; (b) tra- 
jectory in the momentum space; (c) Hikami box along with 
Lyapunov portions of the Cooperon and Diffusons. 

In the quantum problem an amplitude to evolve from 
an initial to a finite point in the angular momentum space 
is given by the sum of amplitudes of all classically allowed 
trajectories passing through these two points. Generi- 
cally, different trajectories come with random and uncor- 
related phases and thus do not produce a systematic in- 
terference contribution. An exception to this rule comes 
from the trajectories having almost (up to k) exactly 
the same geometrical length and thus the same phase. 
This situation takes place if a trajectory contains a self- 
intersection point in the angular momentum space. Then 
another trajectory may exist that is identical to the ini- 
tial one safe for the direction of propagation along the 
loop, see Fig. ^b. The fact that the backward propa- 
gation along the loop is consistent with the equations of 
motion is guaranteed by the time-reversal symmetry of 
the Hamiltonian: 

l^l, 9-^-9, t^-t, H H . (8) 

(compare it with the time-reversal symmetry in the ran- 
dom potential problem: r — > r,p ^ —p,t —t)- It is 
thus easy to see that we are interested in the loops that 
not only have (almost) coinciding initial and finite mo- 
menta: li K Iq, but also (almost) opposite initial and fi- 
nite angles: 9i ~ —9q. The allowed uncertainty is limited 
by the effective Planck constant: (^i — lo){9i + 9q) < k. 

Since such two trajectories have (almost) the same 
phase, they interfere constructively and thus lead to a 
systematic (localization) quantum correction. The prob- 
ability to complete the loop in time t is called Cooperon 
and denoted as C{lo,9o;li,9i;t). With the above men- 
tioned conditions on the initial and finite points one 
may estimate it as C(t) ~ \/Dci/t (this must be mul- 
tiplied by k to take into account the phase area of the 
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allowed uncertainty). This estimate translates (basically 
by the double integration over time) to the —ky/Dd t^^^ 
correctiooiii^ to the classical law: S{P{t)) — 2Dcit . At 
t ^ II = Dci/k^ the correction exceeds the classical re- 
sult and the QKR crosses over to the strong localization 
regime. 

The qualitative reasoning given above repeats identi- 
cally the one employed in the discussion of a particle in 
the field of the "quantum" white-noise random poten- 
tial. However, the kicked rotor dynamics possess a very 
important distinction from that of the white-noise poten- 
tial problem. The latter is the process without a memory. 
Indeed, two classical trajectories that identically retrace 
each other up to a certain point may take completely 
different roots (in particular counter-propagating ones) 
after a single "quantum" scattering event. After complet- 
ing the loop in the opposite directions according to the 
classical random walk dynamics, another single "quan- 
tum" scattering makes the two trajectories to be identi- 
cal again. These two (actually four, since there are two 
trajectories involved) quantum scattering events consti- 
tute the, so-called, Hikami box^^^, denoted by X in Fig. 

[Di. 

Contrary to this scenario, the kicked rotor scattering 
events are purely classical, namely the free rotation of 
the angle. Indeed the trajectory is uniquely defined by 
the standard map, Eq. sequence. If two trajecto- 
ries coincide exactly at some point {On, In), they continue 
to be identical (determined by Eq. (jnj) forever. This 
seemingly precludes any possibility to develop the weak- 
localization scenario, outlined above. The way out of 
this apparent paradox is to recall that the loop may be 
completed not exactly, but rather up to a small uncer- 
tainty: SldO ^ k. This small initial difference is mag- 
nified (more precisely, exponentially increases) upon suc- 
cessive kicking (Lyapunov instability), leading eventually 
to the two counter-propagating diffusive roots. The sit- 
uation is rather similar to the localization in the field of 
the classical large scale random potential (so-called ran- 
dom Lorentz gas). The latter was considered by Aliener 
and Larkin some time agc^^. Due to complexity of the 
Lorentz gas dynamics, they had to introduce quantum 
impurities (essentially a weak white-noise component of 
the scattering potential) to treat the problem analyti- 
cally. The beauty of the QKR is in simplicity of its classi- 
cal dynamics, Eq. ||HJ), that allows to treat the Lyapunov 
regime exactly without involving an artificial quantum 
scattering. 

To proceed in this direction, consider two trajecto- 
ries that initially happen to be at a small distance from 
each other in the phase space: {SO, SI). Taking variation 
of Eq. 0, we find that the angle difference evolves as 
S9n = SOn-i{l + K cos On-i) + Sln-i ■ In the limit > 1 
it leads to SOn = {SO -I- SljK) n^ld (1 + KcosOk) ~ 
{SO + Sl/K)e'^'^, where A is the Lyapunov exponent. For 
if > 1 one finds: A = (ln(ii: cos 6*)) = h).{K/2)^ with 
the residual term ~ 0{\/K). For SOSl « k, the opti- 
mal value of {SO -\- Sl/K) is ^JkjK. It thus takes time 



n = tE = ^ In yjK/k, cf. Eq. (01, to evolve from 
the initial angular uncertainty SO w ^JkjK ^ 1 up to 
SOn ~ 1, when the diffusive motion takes over. Once this 
deviation is reached, the usual diffusive spread of the 
two trajectories takes place. The time-reversal invari- 
ance dictates that aforementioned divergence of the two 
trajectories is preceded by their convergence. The latter 
takes another tE kicks to be completed. The total dura- 
tion of the one-way travel through the Lyapunov region 
is thus 2tE- The entire weak-localization loop construc- 
tion requires two such travels (each including convergence 
and divergence) through the Lyapunov regime. As a re- 
sult, the localization corrections are delayed by (can not 
be developed in time less than) AtE, see Eq. 

Technically there are two equivalent ways to incor- 
porate the Lyapunov region into the weak-localization 
calculations, see Fig. QJ;. One approach, adopted in 
Ref. is to redefine Hikami box to contain AIe scat- 
tering events (kicks) instead of conventional four. Then 
Cooperon (and the Diffusons) is just a conventional diffu- 
sive propagator in the momentum space. In the present 
paper we find convenient to follow the traditional treat- 
ment of Hikami box as consisting of the four kicks. These 
four kicks are treated exactly by multiplying the corre- 
sponding quantum evolution operators. It is the ana- 
lytical expression for Hikami box that dictates allowed 
deviations of SISO product. The Lyapunov regimes are 
delegated to the "legs" of the Cooperon (and Diffu- 
sons). The latter now understood as a propagator that 
includes both Lyapunov-like divergence-convergence of 
the close trajectories along with the normal diffusion 
once a macroscopic deviation between them is reached. 
As we show below, the choice of the four kicks in the 
Lyapunov regime, that coined to be Hikami box is im- 
material. For any such choice the quantum correction, 
linear ink, is delayed by -itE- It is important to mention 
that the time interval < i < 4^^; is protected from the 
higher order loop corrections as well. To demonstrate 
this fact we performed the two-loop weak-localization 
calculation and found the corresponding contribution 
S(P) ~ DciO{t - 6tE){t - QtE)'^/tL to be delayed by GtE- 

The delay is not a bsolutely sharp, but rather is slightly 
smeared by StE = \J\'iIeI^ number of kicks. The 

reason for this smearing is in fluctuations of the exponent 
A. Such fluctuations are due to the fact that one follows 
the Lyapunov instability for a finite number of kicks only. 
(Indeed, unlike classical problems, the minimal deviation 
is limited by k and thus the time to leave the Lyapunov 
regime is finite and may fluctuate between the trajecto- 
ries). Following Ref. QH, we characterize fluctuations of 
the Lyapunov exponent by the other exponent, A2 (for 
QKR with K ^ \ one finds A2 ~ .82). In the case 
Ie ^ 1, cf. Eq. (0), the effect of smearing due to A2 is 
rather small. 

The predicted time-dependent momentum dispersion 
graph is plotted on Fig. [21 The following sections serve 
to quantify the qualitative semiclassical picture outlined 
above. 
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FIG. 2: The momentum dispersion for K = 6.1 and K = 0.6 
- full line; the classical limit ^ 0) - dashed line, standard 
weak localization {tE = 0) - dashed-dotted line; the limit 
A2 0, Eq. - long-dashed line. 



III. CLASSICAL LIMIT: DIFFUSON AND 
COOPERON 



and subsequent averaging over + lo-)/2. Let 

us denote such correlator thereby obtained as 
w), where lj — lj+ — From 
Eq. pU|) one may check that it satisfies 

V(l+,l_;l'^J'_;oj) = e-^^^^ 6i^^i,Ji_j,_ (12) 

+e- {i+\u\i+){i-\u\i'y-T^{i+J'^;i'+J'-;^)- 



The matrix elements {l+\U\l'_l) and {l^\U\l'_L)* may be 
explicitly written as 



{i+\u\i'i)mu\o 



de+ de^ 



cxp 



2^ 



iK COS0+ i0+ , , ' 
k + — 



X exp 



iK cos 6*- 



2k 



(13) 



We proceed now to develop the qualitative considera- 
tions outlined above into an accurate theory of the QKR. 
The essential starting point is the classical diffusion in 
the angular momentum direction of the classical phase 
space: [— tt , +7r] x [—00 , +00] . We'll clarify first how to 
find the classical diffusion from the exact quantum corre- 
lator. One starts from introducing the exact one period 
quantum evolution operator as: 

U = VJ, (9) 
V" exp(iP/2^), J = exp(iXcos^/S) . 

All physical quantities may be expressed in terms of the 
matrix elements of C/", where n stays for the number of 
kicks (time). We are particularly interested in the four- 
point "density-density" correlator, defined as: 

!?(/+, Z_ ;n,n') 

= (/+|J7"e^|;;)(Z_|f7"'e*|;L)*, (10) 

where \l±) denote basis of discrete momentum eigenstates 
oi K — Q quantum Hamiltonian. Note that the forward 
and backward trajectories in this expression are, in gen- 
eral, different. However, it is natural to expect after some 
transient processes, the correlator will be dominated by 
the case of n = n'. This is easy to see by Fourier trans- 
forming the correlator with respect to n,n', passing to 
the frequency ([j+,w_) representation: 

V{1+,1.-J'+,1'_;lu+,lu^) (11) 

n,n'— 



For what follows, it is convenient to introduce the Wigner 
transform representation as: 

1+-I- ^ ' 
xV{1+J_;1'^J'_;lu), (14) 

where we define I = {1+ + /_)/2 and I' = {l\ + l'_)l2. 



A. Frobenius-Perron-Ruelle equation: classical 
kicked rotor 

Let us consider a solution of Eq. (|12|l . As we will see 
below, the integral over 0± in Eq. H13|) are dominated by 
\0j^ — 0- \ ^ k/K <^ 1. This then allows for a perturbative 
expansion of |0+ — 0_ | in the exponent of Eq. (|13|l . In this 
subsection we show that the leading term in such expan- 
sion leads to the classical equation of motion (standard 
map), Eq. ® (so-called "semiclassical approximation"). 
The semiclassical solution thereby obtained is denoted as 
I?o- It should be emphasized that the classical diffusive 
propagator can not be recovered at this stage. To achieve 
this goal, the further approximation must be used, which 
will be clarified in Sec. IIII 51 

The semiclassical treatment employs the following ap- 
proximation for (Z+|C/|Z")(Z_ |J7|Z" )* matrix elements, cf. 
Eq. 1121), 
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{i+\u\ii){i.\u\i'i_r 



27r 27r 

1+ + 1- I'L 



exp 



+ 1 



i 6*4 



2^^ 



2 



(15) 



Let us insert it into Eq. and perform the Wigner 
transform. We also define 9 = {e++9_)/2 and 9' = {9'^ + 
9'_)/2 to simplify the notations. Then with {9+ — 9-)/k, 
{1+ - (/^ and - l"_)/Ti integrated out, 

we obtain: 

Vo{l,9;l',9';uj) = 2nTi 5{l - I') 5{9 - 9' - I) 

+e'^PVoil,9;l',9';u;), (16) 

where P is the Frobenius-Perron-Ruelle (FPR) opera- 
tor, acting on the nearest two arguments from the left 
according to: 

P f{l,9;l',9') 

dh d9i S{l-h-K sin 9i)S{9-9i^l) 

xf{h,9i-l\9'), (17) 

where f {l,9;l' ,9') is an arbitrary function. The 
kernel above implies that the correlator thereby ob- 
tained describes the deterministic motion of classi- 
cal kicked rotor, i.e., standard mapping. Note, that 
in the time representation, Vq is normalized, namely 
(27r^)-i // dld9Vo{l, 9; I', 9'; n) = 1. 

B. Diffusion approximation 

To recover the classical diffusion the procedure 
of deriving FPR equation must be appropriately 
regularizediiiiSi^l. Indeed, in the presence of noises one 
is able to do so and find the proper diffusion constant. 
We shall not follow this procedure here, but rather refer 
a reader to Refs. ITllfT^ 

We shall show below, however, that a regularization is 
consistent with the generalization of Altland's^'* diagram- 
matic method. The latter starts from exact quantum 
density-density correlator, Eq. ifTUI) . and leads to the clas- 
sical diffusion with the correct diffusion constant. The 
basic idea is that for sufficiently large there must be 
some time scale, say Tc, beyond which any classical tra- 
jectory loses memory about its initial conditions. Tech- 
nically, complementing the semiclassical approximation: 
\9n^ — 9n_ I ^ 1 , n = 1, 2, • • • , it is further required that 
1+ = l^ = I and l'^ = l'_ = V at two ends, as well as 
l'_[_ — l"_ = I" for the intermediate variables at multiple 
times of Tc. 

For Tc = 1 the angular memory is lost after every kick 
(cf. the first diagram on the right hand side of Fig.|21b). 
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FIG. 3; Schematic representation for the diffusion approx- 
imation of classical density-density correlator - the solution 
of Frobenius-Perron-Ruelle equation (a), with the self-energy 
(b). The structure of the self-energy is such that the begin- 
ning and ending (after Tc kicks) pairs of the angular momenta 
are the same (denoted by vertical dot-dashed line) , while the 
other pairs not. 



In this approximation Eq. (|12|) is reduced to 
Vo{1,1';lu) 

= Sij,+e'^J2(^\U\n{l\U\iyVo{l'',l'-,u;),m 



where the matrix element {l\U\l"){l\U\l")* is expHcitly 
written as 



{i\u\i"){i\u\i")* ^ {i\j\i"){i\j\iy 

'iK 



d9+ d9^ 



exp 



(cos 6*+ — cos^ 



X exp 



i{9+-9^){l-l") 



d9+ d9^ 



■ exp 



X exp 



27r 27r 

i{9+ -9-){l- I") 



2iK . 9+ + 

sm 

k 2 



■ sm ■ 



(19) 



From this we see (as mentioned above) that the integral 
is dominated by \9+ — 9_\ ^ k/K. Making the change of 
variables: {9+,9^) ^ [{9+ + 9^)/2,9+ -9^ = kip] and 
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integrating out {9+ + 6'_)/2, we simplify Eq. (|19|l as— 



{i\u\i''){i\ur 



dip 
Ztt 



2K kip 

sm — 

k 2 



(20) 

Here J„(x) is the Bessel function of order n. Turn- 
ing to the Fourier representation: 'Df){ip]uj) = 

T.I-V e~"^('~'')l?o(^ - l';^), Eq. (PI gives the classical 
propagator as 

^0(9'; ^) = 1— • (21) 



1 - e'" Jn 



In the limit <C 1 , Kip ^ 1, it is reduced to the usual 
diffuson (Fig.EJ: 

Vo{^;uj)^—^ 2 (22) 

— lUJ + Uci ip 

with the diffusion constant Dci{K) — K"^ /A. 

In Ref. 0, it was shown that the higher order cor- 
relations, namely Tc > 1 [e.g., the second, third, etc. 
diagrams on the right hand side of Fig. ^p] lead to the 
modification of the diffusion constant according to 



Dci{K) 



[1 - 2J2{K) - 2Jl{K) + 2Jl{K) + • • •] 



(23) 

We should emphasize that although the original deriva- 
tion is based on pure classical considerations, it is fully 
compatible with the general formalism developed in this 
paper. For the clarification, we reproduce Eq. (|23|l from 
the exact quantum density-density correlator, Eq. (|12|l . 
in Appendix fXl 

One may wonder whether the dynamics of a classical 
system may be indeed described by the diffusion in the 
long time limit. The answer is known to depend on the 
initial conditions. In particular, it is well-known that 
for any there exist stable islands in the phase space, 
where a trajectory exhibits quasi-periodic motionSiiS. It 
has been estimated, however, that the total area of these 
islands on the phase space is exponentially small in the 
limit K ^ liSiS^. Our approximation, being an expansion 
in powers of 1/K, is bound to lose information about 
these islands. It has also been known that for some 
parametric regions of K, there are some peculiar islands 
(so-called "accelerator mode"^'^^), starting from or near 
which a trajectory will be boosted, faster than the nor- 
mal diffusion's. ( Similar phenomenon exist also in a 
generic modelAS*^. ) Throughout the paper we will stay 
away from these parametric regions, and focus on the 
large K case with the assumption that the initial condi- 
tions do not fall neither into the stable islands nor to the 
accelerator modes. 



C. Classical Cooperon 



1^2-1- — ^2-1 ■ ■ ■ ^ \9n+ — Sn- \ <^ 1. In this part, we turn 
to the discussion of another important limiting case, i.e., 

1^1+ + On-\ ^ 1^2+ + \0n+ + 01-1 « 1 (Cf. 

Fig. El for notations). 



1. Diffusion approximation of the solution of FPR equation 

To find the density-density correlator in this limit let 
us pass to Wigner representation (notice the crucial dif- 
ference with Eq. H14|l ): 



C{1,0;1',9';lu) 
E E expf-|[(?+-r)0-(/V-L)e'] 



(24) 



xCil+J-;l'+,l'_:iu), 



where we define I = {1+ + V_)/2 and V = {l\ + l-)/2. 
The key observation is that the diffuson diagram may 
be retrieved with the bottom (advanced Green function) 
line [Fig. 0] (a)] rotated. This is due to the fact that 
such a procedure simply leads to the time-reversal of the 
series: 6'i_ ^ 6*2- ■ ■ ■ ^ 6'„_ such that it becomes 
—On- —Sn-i- ^ • • • — > ~9i^ . For this reason, we call 
it Cooperon as introduced in Sec.|n] Since |0i+-|-0„_| ^ 
1 , we may proceed along the lines of derivation of Eq. 
from Eqs. (|12|l and H14|) . and arrive at 

Coil,0;l\O';u) = 2TTk5{l-l')5{e -6' -I) 

+e''^PCoil:e;l',e';uj). (25) 

That is, the classical Cooperon is also a solution of FPR 
equation. 



e,+ e,+ 



(a) 



X -x-'-x 
en- 



(b) 



-8i,- 



— 

-e2--ei- 



FIG. 4: The diffusive Cooperon approximation for the clas- 
sical density-density correlator. Typical diagrammatic rep- 
resentation for the classical Cooperon (16^1+ -I- 6n-\,\92+ + 
- <C 1) (a). Rotating the bottom line of (a) re- 
trieves classical diffuson with 9k- replaced by —9k- {k = 
1,2,.- - ,n) (b). 



Starting from the exact quantum density-density cor- 
relator, one may see that the classical diffuson is es- 
sentially related to the limiting case: — ~ 



The only difference with the Diffuson is that instead 
of having poles at kp = 6^ — 6^, the Cooperon has poles 
at kip = 9+ + 9^ (i.e., kip = 9i+ + 9^- = 6*2+ -I- 9{n-i)- = 
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• • • = 6n+ + 6*1-)- One concludes, thus, that the diffusive 
form, Eq. holds in the hmit cj < 1 ,Kip < 1 (cf. 

Fig. 2J for the Cooperon as well. In particular, in the 
case of |0 + 0'| ^ 1 , the averaging over {6 + 9')/2 may 
be performed. Furthermore, if |^ — ^ K , then Cq does 
not depend on the angular momenta. In other words, 

Cq {1,9; I', ~0']uj) = (Co 0; I', -O'-^uj)) f^g^g/y^ 

[Co (w)) , 



_ , (26) 

Ztt —tuj + Ucl f> 

where |/ - l'\ < K, and |0 - e*'] - 1 . 

2. Treatment of the Lyapunov instability regime 

The above general solution for the classical Cooperon, 
Co{l, 9; I', —6'; n) , characterizes the probability for a tra- 
jectory, initiating from {1,0), to end at {I',— 6'). From 
now on we focus on a special case, where SIq = I — 
l';S9o = 0- 9' are such small that \51q\ <^ K, \5eo\ -^l . 
In this part we show that it differs from (Co(a;)) by a 
renormalization factor. 

Without loss of generality, we assume that 
Co{l,9;l',~d';n) evolves from some initial distri- 
bution f {1,9; I',— 6'), bearing the symmetry of 
f {l,9;l',-9') f{l',9';l,^9). Then the formal 

solution of the FPR equation is given by 



Co {l,9;l',^9';n)^ P''f{l,9;l',-9') 



(27) 



To simplify the notations, we rewrite the action of FPR 
operator on Co as: 

P Co {I, 9; I', 9')^ Co {S-^ [I, 9]; I', 9') , (28) 

where S and S^^ are defined as 

S[l,9] = {l + Ksm{9 + l), 9 + 1) , 
S-'^[l,9] = {l-Ksm{9-l),9-l) , (29) 

respectively. Moreover, we introduce the time-reversal of 
FPR operator, Pt, as 



f{l,9;l',9') Pt = jjdhd9if{l,9;h,9^) 

xS{h-l' ~ Ksm9i)S{9i^9' -I) . (30) 

Owing to the symmetry of /, mentioned above, one 
may introduce the following identity: 



Pf{l,9; l',9')^ f{l,9; I' ,9') Pt; 



(31) 



its proof is given in Appendix ^ Applying repeatedly 
this relation to the formal solution, Eq. (|27|l . we obtain: 

Co {l,9;l',~9';n) 

= p"-2"' /(/,6l;/',-6'') } 



= P 



(32) 



for an arbitrary integer number n' such that 2n' < n . 
Consider two nearby trajectories described by {li,9i) 
and {I'l ,9[) , respectively. Their motion is induced by 
S-'^ihO] and 5"[/',-e''] following 



{h,9^) = S-'''[l,9], 
{l[,-9[) ^ S'^[l',~9'] 



(33) 



Associated with the exponential separation of these two 
nearby trajectories, the time, say n^, is defined such that 



\lr,^-l'\^K, \0r, 



1 



(34) 



At such a moment the separation reaches some macro- 
scopic size. After this time the separation experiences the 
usual diffusion in the angular momentum space, while the 
angle difference is uniformly distributed. We substitute 
then n' — into Eq. I|32|l . and arrive at 

Co (/, 9; I', -9'; n) ^ 9 {n ~ 2^^) (Co {n - 2n^)) (35) 

where {Co{n)) is the Fourier transform of (CQ(ti')). Con- 
sequently, the Fourier transform of Eq. (|32|l with respect 
to n is reduced to 



Co{l,9;l',-9';u;) = e' 



dip 
2^ 



Del 



(36) 



In general, ric is determined by the initial condition, 
namely the center of mass: [(/ -I- /') /2,{9 + 9')/2] and the 
initial deviation: {61,69) . One can perform the change 
of variables with respect to general separation {61, 69) ac- 
cording to z = \n\69\, a = 61/69 . In what follows z and 
a are identified as slow and fast variables, correspond- 
ingly. If the initial deviation is small enough such that 
the typical rtc S> 1, the fluctuations of ric at this time 
scale are small. As a result, after averaging over initial 
a, as well as the center of mass, one may cast exppiwric] 
into the renormalization factor (with the logarithmic ac- 
curacy) as (see Appendix lO for details) 



>Vc(2w) = (exp [2icjnc]) 



with 



A 



exp ^iujt 



In- 



c 



2u?X2f 



(37) 



(38) 



in the limit 1^X2/ X'^ <C 1. Consequently, Eq. (p??)]! is re- 
duced to 



Co{l,9;l',-9';u)^kWc{'2uj) 



1 



dip 

271 —ILU + Del P^ 



(39) 
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IV. WEAK DYNAMICAL LOCALIZATION IN 
KICKED ROTOR: ONE-LOOP CORRECTION 

As explained above the weak dynamical localization in- 
volves couplings between the DifFusons and the Cooperon 
[Fig-E]. Therefore one needs a technique to treat two dif- 
ferent kinds of the Wigner transforms introduced for Dif- 
fusons, Eq. 1)14(1 . and for Cooperons, Eq. (|24|l . in a uni- 
fied way. To develop such a technique is the central task 
of this section. We show then that the constructive in- 
terference between two counter-propagating trajectories 
leads to the usual one-loop quantum correction, which is 
a precursor of the dynamical localization. In particular, 
the one-loop correction to the diffusion constant will be 
calculated. 



L I, 



FIG. 5: Sketch of a general diagram, leading to one-loop 
approximation. 



A. Exact interaction vertex 

We will show in Appendix ID II that the one-loop cor- 
rection to the density-density correlator reads as 



St) = 'D (e"^ Vt-i) C (e^'^ V - 1 



-1] V 



with 



C = e"^ VjCo e"^ Vj 



(40) 



(41) 



in the exact quantum operator representation. Here V = 
Vv v. J and Vt = Vj Vv ■ The matrix elements of Vv 
and Vj are explicitly written as 



{1+ M'Pv\l'+ .1'-) = exp 



9+,9.\'Pj\e'^,e'_) - exp 



iK (cos 9+ 



cos 0-) 



x6{6+-6'+) 6{6^-e'_) (42) 

in the representation of the angular momentum and the 
angle, correspondingly. Since we are ultimately inter- 
ested in the long time effects and therefore in the low 



frequencies the e*" factors in Eqs. (|40|l and (|41(l may be 
safely omitted from now on. Thus, we shall not write 
them hereafter. 



1. Minimal wave packet 

In order to calculate ST) explicitly, we consider a gen- 
eral quantity, say Tq = AGB (understanding that A = 
V (Vt - l) , S = [r -ij V for SV), and write it ex- 
plicitly as (for simplicity we omit uj argument) 

I,{l^,l^-J'^,l'_)^J2EJ2E (43) 
in the angular momentum representation, where X is 

x(;3 e^*™^^ h)ll'^ 



§■ cose 


I', 


>" 










,0 


I', 



(44) 



Furthermore, we write down explicitly the matrix ele- 
ments of X, Eq. as 



X 



dOi d9[ (102 d0'2 
where S^in and Sp are defined as 

Skin — 



exp 



{Skin "I" Sp) 



'l ih~l2)+0'Al[~l2) 
'2 ik - h) + 0'^ {I', - I',) 



and 



Sp — K (cos 01 — cos 0'i + cos 02 — cos 02) 



(45) 



(46) 



(47) 



respectively. For the discussions below, we make the fol- 
lowing change of variables: 



<l>={0l-0[)-{02-0'2) 

Sh 

Sh = 



I2 + ^3 ^3 + ^'2 

2 2~ 

^1 + 1'^ 



h^~lm~0[) + i02~0'2)] 



Then Skin may be rewritten as: 



(48) 
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h + I'?, . ^3 + ^2 



2 

h + 0[ 



2 

h + 0'2 



li + I'l ^4 + ^4 



+ 69i 5li + 592 6I2 



[{h - l[) - (/4 - I',)] - ^ [(/2 - I's) ~ ik - 1'2 



(49) 



2 2 

On the other hand, with the semiclassical approximation (i.e. 9i 2 ~ ^1 2) taken into account, S'p may be written as: 



Sri 



-2KS92 sin — cos 



+ 1 



K d> COS 



(50) 



For sufficiently large K, [(6*1 + 9[) - {92 + 6iy]/4 may 
be regarded as a random phase. Therefore, the expo- 
nent becomes self- averaging over this phase. Moreover, 
0. Finally, we employ the conventional hydrody- 



71 2 



namic approximation, i.e. take into account the leading 
term in the if^-expansion only, to arrive at 



exp 



iSr. 



exp 



iSr. 



exp 



iK69i 692 



Jo 



(51) 



Let us insert Eqs. and |(HU, as well as 

Wigner transform of A, B [cf. Eq. itTH) ]. and Co [cf. 
Eq. 12311] into Eq. g3Jl, and integrate out (/i - I'li/h, 
ik-l'i)/^, {l2-l'3)/t:, {h-l'2)/% il'+'l'-)/^ 
and (/)/^. As a result, we find that the semiclassical ap- 
proximation for Xq is 



Ig{l,9;l',9') = V 



A^iW"+^-^,e"+^-^] (52) 



where the vertex operator V is an integral opera- 
tor: V /(/, 9; I', 9'; I", 9"- 5l2,59i) ^ (V/) {I, 9; I', 9') (note 
that the variables: 9", (5/2, and S9i in the function / 
are integrated out) and is defined as 



V/j {l,9;l',9') 
dl"d9" r dSlidS9i 



2Trk 



xCn /" + 



2TTk 



dSl2d692 

592 „_5h 
2 ' 2 ' 



(53) 

X{6h,69i;Sl2,592) 



xf il,9;l',9';l",9";dl2.S9i) 



with 



X{Sh, 59i; SI2, 692) 

^ fK59i592\ 
= Jo T exp 



- {ShS9i + SI2692) 

K 



In Eq. 152() we introduced the following notations 
1 ih + l's , l3 + l'2 



(54) 



I" = 



2 

h+9[ 



2 

h + 



(55) 



while the subscript W stands for Wigner transformation. 

The minimal quantum wave packet posses the uncer- 
tainty: 6liS9i ~ k and 6I2692 ^ k, as indicated in 
Eq. Moreover, S9iS92 ^ k/K. Such quantum 

wave packet determines the initial scale of the deviation 
of two nearby trajectories involved in the Cooperon [cf. 
Eq. 



2. Interaction vertex at the semiclassical level 

Applying the general expression Eq. (|52|l to Eq. H40|l . 
we find 



5V{l,9]l',9') 



(56) 



X P - 1 Po M" - 



2 ' '2 
SI2 



Remarkably, this exact vertex SV is obtained without 
introducing any explicit regularization. In the next sub- 
section we will show that it leads to the weak-localization 
correction to diffusion constant, which is similar to ear- 
lier findings for ballistic quantum dotsiSi In Appendix IeI 
we show that it is possible to introduce some artificial 
quantum disorder to QKR and the results obtained in 
the present work may be reproduced following the for- 
malism of Ref. I23. 
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B. Weak localization correction to diffusion 
constant 

In this part, we show that the exact one-loop quantum 
correction, Eq. (|55|l . may be cast into the conventional 
Hikami box structure with an additional, factor due to 
the Lyapunov region. As a result, the one-loop correc- 
tion affects the diffusion equation through a frequency- 
dependent renormalization of the diffusion coefficient: 
D{uj) = Del + 5D{ijj). That is. 



(57) 



We calculate then the quantum correction 6D{uj) and 
find how it affects the angular momentum dispersion. 



1. Effects of the Lyapunov instability on the interaction 
vertex 

Since our aim is to describe the long-time phenomena, 
we expect that the typical scale of the angular momen- 
tum dispersion, say Lh, is large: Lh ^ K. Indeed the 
angular momentum is randomly spread in the interval 
e [—K,K] in a single kick. It is thus natural to expect 
much broader distribution after many kicks. This con- 
sideration justifies the expansion with respect to K V; 
(hydrodynamic approximation). With such approxima- 
tion, the exact vertex can be cast into the Hikami box. 

With the help of the identity: 'P Vq = Vq "Pt , which 
is proven in Appendix IbI Eq. (|56|l may be rewritten as 



which reflects the time-reversibility and its derivation is 
given in Appendix ^ The remaining procedure is fully 
analogous to calculation of the Cooperon developed in 
Sec. nil CI The only difference is in the boundary con- 
ditions. In fact, the angular deviation of two travelling 
nearby trajectories, 69 reaches a classical size \S9"\ < 1 at 
some point in the phase space, say , 9n^) ■ The later 
evolutions are independent. That is, the two Diffusons 
become self-averaging over the (random) paths connect- 
ing two remote ends, resulting in the factorization of the 
two averaged Diffusons as 



2 ' 



S0i 



( h + —,-9i-—;l, 



x{Vo[h-^-^,-9, + ^-^;l',9' 



(61) 



Taking this boundary condition into account, we obtain: 



51^ 



69" 



;l',9' 



(62) 



dV{l,9;l',9') 



(58) 



SI 



P-l]Vn{ l,9;l" +^-^,9" + ^ 



Since XCq has no dependence on {I", 9"), Eq. H58|l 
contains 



(59) 



To proceed further, we employ the following relation: 
Vo{l,9;l',9')=Vo{l',-9'-J,-9) , (60) 



Note that the two intermediate angular momenta devi- 
ate as 81" jl. Such deviation is unimportant because the 
distribution with respect to the angular momentum fluc- 
tuates over the large scale Lu » X > 51" 12 . In Eq. JHU 
WdI'^) is the same as yVc(a;) except that t'~^ is replaced 
by (with the logarithmic accuracy) : 



In- 



(63) 



where 69 is the initial angular separation of two nearby 
trajectories involved in the Diffuson side of the Lyapunov 
region. It is determined by the minimal quantum wave 
packet, i.e., by X [see Eq. (|^ ] . We then substitute 
Eqs. and into Eq. lf5S|) . and restore the operator 
under the average. As a result, we obtain: 



5V il,9;l',9') 
dhd9i 



2Trk 



Vn I 



SI" 



(64) 



Ft - 1 



P-l]Vo{h-^-^,-9, + ^-^;l',9' 
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where 



V 



dShdSOi f dSl2dS92 



2Trk 



2Trk 



WB{2u;)X{Sh,Sei; 512,662) 



Sh 

T 



S9_2 
2 



(65) 



Here V may be regarded as the renormahzed inter- 
action strength. To further calculate it, we substitute 
Eq. (jSnj into Eq. (|^ . As a result, V is found to be 



V = T{lu) 



dip 



where 

r(c.) 



27r —iuj + Dcif^ 
dShdSOi f d6l2dS92 



(66) 



2TTk 

dShdSOi 
2'Kk 



2nk 
d5l2dS02 



Wc (2t^) Wd i^Lu) X 



2Trk 



X Jo ( ) exp 



-{5h59i+5l2592) 



X exp 



2w2A2 
2iuj ) \zi + Z2I 



(67) 



with 



zi = In ^59l+5ll/K'^ , Z2 = In ^59l + 5lllK' 
Rescaling 59x^2 and 51x^2 as 



59x,2/VWk, 5h^2^ Sh,2/VkK 



0Cl,2 

leads to 

r(cj) = exp ( AiujtE — 
where F{uj) is 
F{u;) = 



(68) 
(69) 



4:Uj'^\2tE 
A2 



F(c.) , (70) 



dlid9i f dl2d9 



2-K J 2tt 
xexp[i {Ii9i+l292)] 



Jo i9i92) 



X exp 



2iLj — 



\zi + Z2I 



A2 



Ji = \nJ9i+li, Z2=\nd9i+l 



(71) 



Here Ie = {t^ + /2, which is Eq. (O. Since zi,2 1, 
in the limit lo, lo\2/X^ ^ 1, the last exponent in Eq. H71|l 
may be considered to be 1. As a result, F{ijj) — 1. Thus, 
we obtain: 



T{uj) — exp AiuotE — 



4cj2A2ti; 

A2 



(72) 



We point out that the position of minimal wave packet- 
Hikami box-can not be exactly located within the Lya- 
punov region. Indeed, this is reflected in the fact that 
the total duration within the Lyapunov region, i.e., AIe 
actually does not depend on the exact boundary between 
Cooperon and Diffuson. Such feature originates from 
the chaotic nature of the classical motion in the phase 
space. At each full travel, the initial deviation, 5li,2 
{59i 2) with respect to the reference trajectory expands 
in the backward (forward) time direction. Eventually 
5li 59i {5I2 592) reaches some classical action K (the typ- 
ical scale of the classical action) . Therefore, the total 
duration is 



4t, 



A 5h59i 



-In-^ 

A ^ 612592 



(73) 



Taking into account the uncertainty relation: 5li59i k, 
6I2592 ~ k , w e find the total duration to be AtE — 
4:\-Hn./Kjk. 



2. Frequency-dependent diffusion coefficient 

The renormalized interaction vertex, Eq. H64() may be 
further cast into the conventional (diffusive) Hikami box 
upon further simplifications. As discussed above, the 
hydrodynamic expansion may be performed because of 
K/Lh <C 1. This allows the further simplification of 
Eq. (|64|1 . One finds to the first order in the hydrody- 
namic expansion 
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51" 

= {Vo{h- — + Ksm 



69" 51" 



„ 59" , 51" 
h - ^ + li - — 



59" 51" 
l + /fsin(0i-^ + /i-^) Vz, 







-Do 















2 2 

(561" , (5/" 
1 - ^ + ^1 - ^ 



59" , (5?"' 



(74) 



The last line results from the fact that (I?o) has weaker dependence on the angle compared to the sinusoidal term for 
sufficiently large K. Similarly, 



51" 

Vo [ i,e-h + —,Oi 



\ -\- K sin 



59" 51" \ 



51" 59" 51" 
Vo{l,9;h + —,9, + — + h- — 



(75) 



On the other hand, by shifting the overall angle factor, 



dlid9i 
27rS 

dlid9i 
27rS 



Vo I 



51^ 



+ ^1 



51" 



Vo[h + ^-^,-9, + ^-^-,l',9' 



51" 

Vo[l,9;h + —,9, 



51" 59" 
Vo ( l,9;h + —, + — 



51^ 



(76) 



This arises from the uniform distribution with respect to the common angle in the product of the two averaged 
Diffusons. 

We then substitute these two expansions: Eqs. (|71|l and (|75|l . as well as Eq. (17^ into Eq. For sufficiently 

large K, the sinusoidal term is quasi-random and may be averaged over the angular region [0,27r]. As a result, the 
linear term in the hydrodynamic expansion does not survive upon this averaging and the second order term must be 
kept. Finally we obtain: 



5V (/, 9; 9') ^ V I /sin (^i - ^ I sm i 



51" 

Vo[l9-l[ + —,9^ 



Vo I'l 



51" 



'i 'i 
„ 59" 



(77) 



r 



Since 59" < 1 - 6ii, 



( sin 



59" 



sm Wi 



59" 



(sin^ 9i) - 2D,. 



(78) 



On the other hand, the Diffuson is smooth over the scale 
~ K, hence we may also drop out 51" / 2 in Eq. (|77|l . 
Finally 5V is cast into 

5Vil,9;l',9')^2V j'^^D^iVvVv; (79) 

X [{Vo {I, 9; l{,9,)) {Vo {I'l, ~9^; I', 9'))] \v^^v(=i, • 

The diffusion coefficient appearing in Eq. H79|) is Dd — 
We argue that if more kicks are reserved for the 
Hikami box, i.e., {P — l)",n > 1, then the diffusion 



coefficient will acquire the same higher order corrections 
as Eq. 

With the angle averaged out diffusion is retrieved. As 
a result, the one-loop correction, Eq. ifT^ . is simphfied 

as 



SVo{l,l') ^k-^ / dhVD. 



V 



x[Vo {lJ[)Vo {l'l,l')]\v-v^=i. 



(80) 



Denote the Fourier transform of Vo {I, I'), with respect to 
l-V ds Voi(p;uj). Substituting Eq. (jHU into Eq. (jSOl) we 
find the Fourier transform one-loop correction: 



5Vo{ip;uj) 



kT{iu)D,i 
{-iu! + Del •^'^f 



1 



TT —iuj + Del 0^ 



(81) 
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which leads to the one-loop quantum correction to the 
diffusion coefhcient as 



SD{u;) = -^T{lu) 



(82) 



3. Dispersion function 

One may express the time evolution of the angular mo- 
mentum dispersion as: 6{P{t)) = {{l{t) — 1{0))^) in terms 
of the frequency-dependent diffusion coefficient. In fact, 
by averaging over the angle, we may write S{P{t)) as 



d(p 



(83) 



Substituting Eq. H57() into it, we obtain: 



duj 1 



Diio). 



(84) 



For sufficiently large K one may ignore the fluctuation 
of A, i.e. put A2 = 0. Consequently, in the leading order 
in k the momentum dispersion is found to be: 

6{l\t)) = 2D,it - e{t - AtE) {t - AtEf'^ , (85) 

where 6{t) is the step function (long-dashed line on 
Fig. 121). The singularity at t = 4^^; is rounded by the 
Ehrenfest time fluctuations-arising from finite A2 (full 
line on Fig. 01. Substituting Eqs. ((TH and into 
Eq. H84|) . we arrive at 



6TT 



(86) 



3/2 



e 16 



5tE J V <5*-E 

where SIe = \/ ^2iE / ■ As a result. 



37r/64 



5tf 



where /(O) = 1 and 
8V2r(i) 



(87) 



fix) = { i^(t) 



5/2 ^-x^'/m for X > 1 



for — a; 1 . 



This result completes the calculations of the one-loop 
weak-localization correction. 



V. WEAK DYNAMICAL LOCALIZATION OF 
QKR WITH BROKEN TIME-REVERSAL 
SYMMETRY 



To exploit farther similarities and differences of the 
dynamical localization of the QKR and the Anderson lo- 
calization we discuss here effects of breaking the time- 
reversal symmetry (TRS). In case of Anderson localiza- 
tion in a random potential the TRS may be broken, for 
example, by a static magnetic field. The latter provides 
different phases for clock-wise and anti-clock-wise prop- 
agating trajectories, destroying thus the systematic in- 
terference correction discussed in Section II. It does not 
ruin, though, the Anderson localization completely. In- 
deed, higher order corrections (the minimum possible is 
the two-loop one) may be interpreted as interference of 
trajectories travelling the loops in the same direction only 
(Diffuson only diagrams with no Cooperons). The static 
magnetic field does not affect such diagrams. As a result, 
the Anderson localization exists even in this case, albeit 
with somewhat larger localization length. 

It has been thus of long interest to show that an anal- 
ogous phenomena exists for the QKR as well22i2li^2ii2iii. 
Our additional motivation comes from consideration of 
the Lyapunov regime and its sensitivity to the TRS 
breaking. In particular, the one-loop (TRS invariant) 
correction was found to be delayed by A^e- Does this 
time interval remains to be protected against perturba- 
tive corrections in higher loop processes ? Is the delay 
time still the same ? These questions are of particular 
interest if and when the leading one-loop correction is 
destroyed by TRS breaking. 

To answer these questions we investigate the model, 
described by the following Hamiltonian 



cos 6* 5{t-2n) 



-hcos {e + ^^ 5{t- (2n+ 1)) 



(89) 



Here the time-reversal symmetry is broken for generic 
$ except for $ — 0, tt . First, we analyze suppression 
of the Cooperon (and thus the one-loop diagram) aris- 
ing from We then calculate the two- loop correction. 
In contrast to the one-loop correction, the two-loop one 
is robust against <& because it contains (among others) 
a diagram without the Cooperons. If the TRS is bro- 
ken, the long-time correction is given by the two-loop 
diagram depicted in Fig. E?^''*^''*'*. It differs from Fig. Q] 
in that the two interfering paths propagate together in 
the same direction, except inside the Hikami box, where 
they switch from one to the other. Due to such geometry, 
the two-loop correction is not sensitive to <& . It requires 
three successive travelling through the Lyapunov region 
(each taken 2tE time). We show thus that the weak lo- 
calization correction, given by this diagram, is delayed 
by Gifiii. 
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FIG. 6: The leading quantum correction to the density- 
density correlator in the absence of the time-reversal sym- 
metry: (a) two-loop weak localization diagram; (b) its image 
in the momentum space; (c) Hikami box. 



A. Suppression of one-loop correction 

In the modified KR model, Eq. (|89|l . the period is dou- 
bled and includes two kicking. The effective kicking op- 
erator appearing in the one-step evolution operator U [cf. 
Eq. l(T^ ] is replaced by 











exp 


I— cosd 

K 


exp 


2k 



X exp 



K 
i— cos 

K 



(90) 



In the angular momentum representation, the matrix el- 
ements read as 



E 



27r J 27r 

xe^e-^«i('^-'+)e^^°^(^i+*', (91) 
while its complex conjugation is 

xe-^ e^''i(''i-''-)e-^^°<^i+*) . (92) 

First we investigate the effects of $ on the diffusive 
parts of Diffusons and Cooperons. To simplify the discus- 
sion we assume that K is sufficiently large, implying that 
the memory about the angle is lost after a single kick. 
To find the self-energy of the Diffuson, we replace the 
kicking operator J in Eq. (|19|l with V'j . Then we insert 
Eqs. H91|) and (|92|l into it, putting li — l[ . Consequently, 
the diffusive pole is retrieved as 9^ — 6^ — 6i — 9[ = kip , 



and the self-energy of the Diffuson is 

I'l-pjini'i-pji'")' 

This expression demonstrates the period doubling (Bessel 
function is squared) and implies that the Diffuson is not 
affected by $ . 

For the diffusive Cooperon, the self-energy is 
{l\Vj\l"){l"\'Pj\l)* ■ Inserting Eqs. ^ and ^ into it, 
and putting li = , we find that the Cooperon has dif- 
fusive pole at 



and the self-energy is 

{i\'p'j\n{i"\'p'j\i) 

kip 



X Jo 



2K 

— 2 



^0 



2K 



kp, 

2tt 

kp 

~2 



(94) 



(95) 



In the derivation above, we used the fact that for suffi- 
ciently large K , {6+ - 6'_)/2 and (6*1 - e[)/2 is quasi- 
random and the self-averaging may be performed. With 
the Fourier transform with respect to time, the self- 
energy leads to the diffusive Cooperon of the form: 



(Co(/,0;/', 
dp 



(96) 



27r 



p^ 



{p + 2^/kf 



in the limit: Kp, K^/k < 1 . Here \l - l'\ < K,\e - 
0'\ ^ 1 . According to Eq. I|9H|I the relaxation time of 

the Cooperon is r$ = (Dc/^'V^^)"^ . At $ > k/y/D^ ^ 
k/K , T$ ^ 1 the Cooperon is completely suppressed. 

The propagation in the Lyapunov region involves de- 
terministic motion, which, is also affected by $ . There- 
fore the Lyapunov exponent A, as well as its fluctuations 
A2 acquire some ^-dependence. However, since the prop- 
agator has no dependence on the center of mass, the func- 
tional form of yVc,D(2ti^) remains unchanged. Thus, the 
procedure of Sec. lIVBl can be employed to show that the 
dispersion function is given by: 



Ak^jD, 



cl 



d{l'{t)) = 2D,it 



where the function h{x) is 



t-itE 



h{x) 



dyi 



dy2 e 



(97) 



(98) 



Notice that here is a function of<I>, i.e. i£;(<&) = 
A(<i>)~^ In ^jK/k . For simplicity, we neglected the fluc- 
tuations of A(<i>), i.e. we put A2 = 0. In the region 
t — AIe 3> T$ , the one-loop correction is exponentially 
suppressed. One is required, thus, to consider the higher- 
loop corrections. The two-loop correction (Fig. gives 
the leading weak dynamical localization correction. 
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B. Two-loop correction 

In principle the technique developed at the one-loop 
level can be employed to treat the two-loop case. How- 
ever, this is technically quite involved and is not discussed 
here. To read out the frequency-dependent diffusion coef- 
ficient in an economical way, let us renormalizc the stan- 
dard results of the weak localization^ with an appropri- 
ate t_E-dependent factor [this procedure indeed is trans- 
parent at the one-loop level, Eq. I|82|l ] . The renormaliza- 
tion factor for the two-loop geometries (see Fig. were 
calculated in Ref. in the context of the Lorentz gas. 
Adopting the analogy between the Lorentz gas model and 
the QKR, verified above on the one-loop level, one finds 
for the two-loop frequency-dependent correction to the 
diffusion coefficient 



5D{uj) = -2k^Dcir3{uj) 
whereSSiii 



dip 



2tt {—iuj + DciLp^) 



Tj,{ijj) = exp ( QiuutE 



9uj^X2tE 

A2 



, (99) 



(100) 



As a result, the leading correction to the momentum dis- 
persion in the case of broken TRS is given by: 

Again we ignore A2 for simplicity. To develop the two- 
loop geometry, a minimal quantum wave packet must 
take time t^; to expand into some macroscopic size, and 
vice versa. Within the logarithmic accuracy appear- 
ing here is the same as that in the one-loop correction. 
Therefore, the duration for a full travel through the Lya- 
punov region remains the same as the one-loop case, 
namely 2tE time. The two-loop geometry involves three 
successive visits. In the diagrammatical language, each 
leg of the 6-leg Hikami box (see Fig. Elb) lasts time i^. 
Thus ,the weak localization correction, given by this di- 
agram, is delayed by 6tE—- 



VI. OBSERVATIONS OF 
CLASSICAL-TO-QUANTUM CROSSOVER IN 
REALISTIC DRIVEN SYSTEMS 

In this section, we discuss some possibilities for exper- 
imental observations of the predicted ig-dependence of 
the classical-to-quantum crossover. The quantity to be 
measured is the dispersion function 5{P{t)). 



A. Energy growth in cold atomic gases 

In the 90's unprecedented degree of control reached in 
experiments with ultra-cold atomic gases*'^ allowed to in- 
vestigate various fundamental quantum phenomena. The 



advent of laser cooled atomic gases and standing wave op- 
tical pulse a^i^i^ has opened the door to study quantum 
chaos experimentally. In an insightful paper^ Graham, 
Schlautmann, and ZoUer pointed out that atom optics 
may serve as a testing ground for quantum chaos. Shortly 
later, the idea came into realization with sodium atoms 
being cooled and trapped using the magneto-optical trap, 
subjected to a phase-modulated standing wave^^. Later 
on, a realization of the QKR in atom optics was ac- 
complished with the phase-modulated standing wave re- 
placed by a pulsed standing wave^. 

In an atom-optical experiment, typically 10^ sodium 
or cesium atoms are trapped and cooled down to 10/iK 
using the conventional magneto-optical trap. After 
turning off the trapping fields, two linearly polarized, 
counter-propagated optical beams with the frequency ujl 
are switched on, creating a spatially periodic potential: 
Vb cos{2kLx) . Here — uj^/c is the laser wave number. 
Such optical lattice is controlled by the acousto-optical 
modulator as a pulse sequence with a profile f{t) . The 
pulse length Tp may be much smaller than the period 
T . The atomic cloud, exposed to this pulsed optical lat- 
tice thereby, experiences a series of kicks. The evolution 
of the atomic momenta distribution is monitored after a 
certain number of kicks. 

In experiments, the laser detuning = lol — ujq 
from the resonance frequency luq is large compared to the 
excitcd-statc decay rate. The dipole force due to Stark 
effect leads to the spatially-dependent shift of atomic lev- 
els. This results in an effective periodic potential imposed 
on the atomic cloud. One may model the center-of-mass 
motion of atoms with the single-particle time-dependent 
Hamiltonian as 



n2 



H 



P 
2m 



N 



Vo COS {2kLx)Y, fit -nT), 



(102) 



ji=0 



where m is the atomic mass, and fci is the laser wave 
number. The effective potential, Vq is determined by the 
maximum Rabi frequency fl as Vq = hfl^ /SA^. The po- 
sition X, and the momentum p operators are canonically 
conjugated: [x,p\ = ih. The effective kicking strength 
K, and the Planck constant k may be expressed then as: 



K = 8uJrT'^TpVo/h , ^ = 8 UJrT 



(103) 



with the recoil frequency cUr — hk\/2m . With the rescal- 
ing: X e = 2kLX, p ^ I = {k/2hkL)p, t ^ t/T, 
fit) fit)/Tp, H -> {kT/h)H , one casts the Hamilto- 
nian into Eq. (P) in the limit Tp/T (with K fixed)''&. 

So far many experimental efforts^ have been focused 
on the parameter range where k ^ 2, and thus is « 1. 
The dynamical localization has been observed as the sat- 
uration of the time-dependent momentum distribution 
width (i.e. energy absorption). To extract accurately 
the i^-dependent crossover one needs a large separa- 
tion between the relevant time scales: 1 < Ie < . 
This requires to decrease k down to 0.1 — 1 , which we 
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hope to be soon within the reach for cold atomic gases 
experiments2Si2^. 

B. Charge fluctuations on Josephson grains 

The experimental realization of the QKR may be 
also feasible in experiments involving small nonequilib- 
rium superconducting systems. One example is a small 
superconductive dot in contact with a bulk supercon- 
ductor through two Josephson junctionsiS. The bare 
Josephson coupling, Ejq, is modulated via the external 
magnetic flux threading the SQUID loop^°i^^i^^: Ej = 
Ejo\coa{nBAioop/^o)\, where Aioop is the area of the 
SQUID loop, and <I>o = h/2e is the superconducting flux 
quantum. If B is modulated in a meander way, with the 
pulse length much smaller than the period T, the system 
may be modelled by the Hamiltonian: 

^ (Q-cvy _ 

H = Ejcose 22^5(1 ~nT) . (104) 

n 

Here 9 , and Q are the relative phase of the superconduct- 
ing order parameter on the grain, and its charge, corre- 
spondingly. They are canonically conjugated: [0, Q] = 
2e i. In Eq. H104|l C is the capacitance, Vg is the gate 
voltage, and Ej is the time-average Josephson coupling. 
Making change of the variables: Q ^ I = hQ/2e 
and rescaling the relevant quantities as t — > t/T ,H — > 
SEcT^H/ti^, and CVg Vg = hCVg/2e [E^ = eV2C), 
we cast the Hamiltonian above into QKR: 

H = ]^{l~vg^^ ~Kcose^6{t-n). (105) 

n 

Note that the sign difference in the kicking term is imma- 
terial. Here the effective Planck's constant and kicking 
strength are 

n = 8EcT/h , K = 8EcEjT^/h^ , (106) 

respectively. 

The charge fluctuations are described by the charge 
dispersion: ~ 5{P{t)) . For sufficiently large K and 
t < 4t£; it is expected to increase linearly in time. At 
t ^ ^tE, it should deviate from the linearitjiSS. At longer 
time, t ^ AtE the t'^/^ power-law correction develops fol- 
lowing the conventional weak localization theorjiiiiSiiS. 
This signals the onset of the localization phenomena. 
Eventually at i ^ Ddj'k'^ ^ ts the charge fluctuations 
saturate and do not grow any more upon further kicking. 

C. Charge fluctuations in superconducting 

nanocircuits 

Recent work24 suggested another kind of time- 
modulated small superconducting system. It is proposed 




FIG. 7: The scheme of the driven Josephson grain: a su- 
perconducting dot (left) is coupled to a bulk superconductor 
(right) via a SQUID loop (middle). The flux $ piercing the 
SQUID loop is time-dependent. It effectively modifies the 
bare Josephson coupling Ejq of the two junctions. 



that a mechanically driven superconducting single elec- 
tron transistor (SSET) may serve as a realization of the 
QKR. The system is based on a Cooper pair shuttle - a 
small superconducting island, periodically travelling be- 
tween two macroscopic superconducting leads with the 
phases 0l and respectively. Twice during each pe- 
riod, T , the shuttle meets one of the leads, experiencing a 
sudden Josephson coupling. The average coupling energy 
Ej is assumed to be much larger than the charging energy 
Ec- If the two leads are far enough from each other, the 
island never couples to both leads simultaneously. If the 
switching time is short, the time-dependent Josephson 
coupling may be mimicked by delta-pulses. Therefore, 
one may model the system with the following Hamilto- 
nian: 

H = - Ej [cos eT5{t~ 2nT) 

n 

+ cos(6l-H$) r(5(t- (2n+ 1)T)]. (107) 

Here is the relative phase of the superconducting island 
with respect to the right lead. Remarkably, the phase 
difference across the two superconducting leads, $ = 
4>R — (f'L breaks the "time- reversal" symmetry, Eq. 
The effective Planck constant and kicking strength re- 
main the same as in Eq. Ij 106(1 . With the same rescaling 
as in Sec. IVI Bl Eq. (|107|l is rewritten as Eq. 

The classical-to-quantum crossover is reflected in the 
nonequilibrium charge fluctuations of the superconduct- 
ing island. In the case of $ = 0, the situation is the 
same as Sec. IVIBI In the presence of a small phase bias 
across the two superconducting leads, i.e., |<I>| ^ ^/K , 
the charge fluctuations are described by Eq. (|97|l . For 
larger phase bias: \^\ > ^/K , the charge fluctuations 
are given by Eq. (|101|) . 
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(a) 




(b) 



T 2T 3T 4T 5T 6T 7T 8T 

FIG. 8: A superconducting shuttle periodically travels be- 
tween two superconducting leads with the phase 0l and <j>R , 
respectively (a). At every other period 2T , the shuttle ex- 
periences a sudden Josephson coupling with the left (right) 
lead. 



VII. WEAK DYNAMICAL LOCALIZATION IN 
PRESENCE OF NOISE 

Due to the quantum interference nature of the weak 
dynamical locaUzation, the effect may be strongly sensi- 
tive to noise. Indeed, noise effects are known to be of im- 
portance for the dynamical localization^Si^i^^iS. So far 
both experimental and theoretical studies have been pri- 
marily focused on < 1. Below we shall consider how 
an external noise affects classical-to-quantum crossover 
at time scales t ^ Ie- To this end, we will employ the 
technique of Sec. llVl to investigate a noise-sensitivity of 
the weak dynamical localization Cooperon and Hikami 
box. 



A. Phase noises 

In proposed experiments involving small superconduct- 
ing systems, the phase may fluctuate due to the influence 
of the dissipative measurement circuit. For the analytical 
treatment we focus on the particular kind of noises - the 
Gaussian phase noise. That is, the phase is assumed to 
fluctuate randomly in time: — > + Cn ■ Here C„ = Ct 
is the noise that is assumed to be uncorrelated at differ- 
ent kicks (CnCn') = for n 7^ n'. At a given kick the 
random phase C is supposed to be drawn from some pe- 
riodic distribution function P(C) = (27i')~^ X^m -^"i^""^- 
For simplicity wc shall assume that 

P„ = e-'"^™' , (108) 
where a^p characterizes the strength of the noise. 



1. Strong phase noise 

The strong noise limit is characterized hy — > oo, 
and thus P(C) = l/27r being uniform distribution on 
C G [0, 27r] . We show first that the classical diffusion 
(without localization) is restored in this limit. The quan- 
tum density-density correlator satisfies (cf. Eq. (|10|l ): 



V^{l+,l_;l\,l'_) 



(109) 



^ {i+\u\il){iAW^)*vai,i--.i'+,i'^ 



where the long bar stands for the average over the phase 
noise with respect to the uniform distribution, i.e.. 



{i+\u\il){i.\u\l'L) 
d( de+ de^ 
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(110) 



X exp 
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2k 



iK 

ie+ 



2n 



cos {d+ + C) - COS (6l_ + C)] 

f {1+ - I'D 5i+-i-..vi-i'i. ■ 



In the last line /(/) is defined as 

/(0 = 



dtp (2K kLp 



(111) 



With the definition: I = [1+ + /_)/2, V = [V^ + l'_)/2, 
and Al = 1+ — ^ the solution of Eq. H109|) may be 
formally written as 



'ilAl 



(112) 



exp 



2k 



'JLl' 



E 



E 

n— 1 h,--- ,ln-i 

Xf{l-h)f{h~h)---f{ln-l-l') , 



exp 



iAl^, 



k=l J 



where Z„ = I' . As a result of the averaging over noises, 
V,^ has no C-dependence. In the large K limit, we expect 
the solution to be independent of {l + l')/2. Thus, taking 
the average over {I + l')/2, one finds: 

oo 

mu')) = ^M'+E E (113) 
^f{i-h)nh~i2)---nin-i-i') ■ 

Passing to Fourier representation, we find that it is 
nothing, but the classical diffuson, Eq. H21|l . It is 
worth mentioning that in the strong noise limit the dif- 
fusion constant is exactly K'^/A. This is due to the 
fact that higher order corrections:^^ as well as quantum 
renormalizationtSi^ result from long-time correlation ef- 
fects (cf. Appendix . The latter is completely de- 
stroyed by the strong noise. 
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2. Weak phase noise limit 

We turn now to a more interesting case, relevant to 
the context discussed here, i.e., ti^ ^ 1, where the noise 
only slightly suppresses the weak dynamical localization. 
For (^^) ~ (T<^ <C 1 , cos{(p + C) ~ cos — C sin ip. 

Let us concentrate on the effect of the noise on the 
Cooperon. Upon every kicking, Cooperon acquires an 
additional phase: K{(t cos 0+ — cos 9^)/k . Recall that 
t , t' are counted from the two opposite ends of the loop 
trajectory, and t+t' = Tt with Tt being the total duration 
of the loop. Since C,t and C,t' stand for noises at different 
moments, they are uncorrelated. Averaging the phase 
factor over them leads to the exponential suppression of 
of every step of the Cooperon ladder: 



As a result, the intermediate intensity noise {k/KY ^ 
cr,p <C 1 acts, to much extent, as a TRS breaking pertur- 
bation. It suppresses the Cooperon corrections, leaving 
the Diffuson ones (the simplest being the two loop one 
Fig. intact. 

Above we find that the Diffuson-Cooperon coupling- 
minimal wave packet, is suppressed by large enough 
phase noises. This picture is naturally expected to be 
applicable for higher order interaction vertex also. The 
latter is essentially responsible for the onset of the weak 
dynamical localization in Diffuson-only diagrams. This 
picture may be considered to be the precursor of the 
restoration of diffusion for larger noises. 



exp I — - — cos 6^4 



exp 



iKCf 



COSW- 



4^2 



(114) 



(/o (x) is the modified Bessel function) . This implies that 
the Cooperon is suppressed as 



Co(t) ->Co(t)e 



(115) 



with the dephasing time defined as = 2^"^ / [K'^a^] . 
For cr^ > {%/KY , ^ 1 , the Cooperon mode is sup- 
pressed completely. It does not mean, however, that the 
classical diffusion is restored for such small a^. Indeed, 
the higher order loop corrections, that include the Diffu- 
sons only (cf. Fig. may still survive such level of the 
noise and lead to the dynamical localization. To verify if 
this is the case, one needs to study the effect of the week 
noise on Hikami box. 

Starting from Eqs. (|48|) and H55|) . we observe that S6i 
and 9" are independently shifted by the noise, i.e., S6i 
S0[ = S01 + Ci , 0" 9" + C,2- As a resuh, instead of 
Eq. if^H) . X is given by 



X{5h, 59i- Sh, 502) 



(116) 



sin 61" X{5h, S0[; Sh, 692) . 



exp 



Upon averaging over C, Eq. (|116|l leads to the the expo- 
nential suppression of the minimal quantum wave packet 
(Hikami box) as 



exp 



2k 



XiSh, 60[; SI2, 602) . (117) 



From here we see that the weak phase noises {cr^ <C 
1) does not substantially affect Hikami box. Indeed, 
the interaction vertex is significantly suppressed only at 
CTip ^ 1 . This means that the effect of phase noises on 
Hikami box may be ignored compared to the dephasing 
of Cooperon. 



B. Amplitude noise 



In cold atoms experiments, the optical pulse power 
may fluctuate with time and thus the series is not per- 
fectly periodic. This then leads to the noise in the kick- 
ing amplitude, i.e. the stochastic parameter K is re- 
placed by K -\- r]n, where r]n — rjt is a random ampli- 
tude fluctuation. The effects on dynamical localization 
of such kind of noise have been under intensified exper- 
imental investigations^Si^^iS. A central issue addressed 
is, whether the dynamical localization is destroyed com- 
pletely by the noise. We concentrate here on a weak noise 
limit. 

To simplify analytical estimations below let us as- 
sume that r]n is the white-noise, Gaussian noise, namely 
iVnVn') — cfK^nn'- Here (Tk characterizes the strength of 
the noise. Below we consider the limiting case of ctk ^ 1, 
where the amplitude noise only slightly suppresses the 
weak dynamical localization. 

For the Diffuson the one-step quantum propagator 
acquires an additional phase as ryt(cos0+ — cos6'_). 
Here 0+ , 9- stand for the phases of retarded/advanced 
Green's functions at the kicks. Passing to the semiclassi- 
cal limit, i.e., \9+ — 0-\ <C 1 and averaging over the noise, 
we find that the Diffuson is affected via renormalization 
of Del by a small correction 6Dci ^ ctk <C Dd (A de- 
tailed study in this direction was presented in Ref. |^). 
This correction is not responsible for the destruction of 
the dynamical localization. It is the effect of the noise on 
the Cooperon and Hikami box that eventually may lead 
to the restoration of the classical diffusion. 

In the case of the amplitude noise the effect on the 
Cooperon is fully analogous to the case of the phase 
noise, i.e. dephasing. Upon every kicking, the Cooperon 
acquires an additional phase: (774 cos 0+ — rif cos 9 ^)/k 
with t , t' counted from the opposite ends respectively, 
such that t + t' = Tt with Tt being the total duration of 
the loop. Upon averaging over the independent noises rjt 
and rjt' , this phase factor leads to the exponential sup- 
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pression of a single step of the Cooperon as: 



exp [ — cos( 



exp 



exp 



2^2 



^2 



(118) 



This impHes that the Cooperon is dephased as Eq. pi5|) 
with the dephasing time being tk = /(Jk ■ For 
^ , tk ~ 1 , the Cooperon mode is suppressed 
completely thus only higher order terms (Diffusons only) 
may be responsible for the dynamical localization. It re- 
mains open to estimate the noise amplitude that destroys 
Diffuson-only weak dynamical localization. 



C. Effects of finite depiiasing time 

We saw above that every realization of the QKR may 
involve various noises in realistic experimental environ- 
ments. As a result, there exist various dephasing mech- 
anisms. The effective dephasing rate is the sum over all 
dephasing rates, namely 



(119) 



To simplify qualitative discussions, in this part we 
focus on effects arising from weak dephsing such that 
He ^ T^- First we show that the effective Ehrenfest 
time is shortened. In fact, the additional weak noises en- 
hance the rate of angular deviation spread according to: 



— (5r(t) = 2A(5r(t) , 
at T, 



1 



(120) 



Recall that the first term on the right hand side results 
from the Lyapunov instability. The solution is easily 
found, i.e., Se'^it) = [S0i,2 + (2Ar^)-i]e^* - {2Xt^)-\ 
Here S9i,2 are the initial angular deviation of Diffu- 
son/Cooperon, respectively. At t^* (or t'-^*), 69 ^ 1. 
The solution then gives 



t 



D,C* 



In 



1 



^Sel^ + (2Ar^)-i 



(121) 



Thus, the effective duration of a full travel through the 
Lyapunov region is found to be 



< 



In 



^J{sel 662)^ + {sef + sej) (2At^)-i + (2At^)- 
1 



In 



In- 



\S9i ,5021 + (2At^)-i 
1 



K 



(2At^ 



= 2n 



(122) 



The effective Ehrenfest time t*^ <tE- 



With t*^ substi- 



tuted into the renormalization factors Wc and Wd , and 
taking into account the dephasing of diffusive Cooperon, 
the one-loop correction, Eq. I|82|) . is modified as 



5D{uj) = - 



Wd i2u) Wc i2uj) 
dip 



(123) 



Taking the Fourier transform with respect to ui, one find 
the momentum dispersion to be 



-= 2Dr,t- 



(124) 



Again A2 = 0. According to Eq. (|124|l the quantum cor- 
rection is suppressed at t > -I- . 

One should keep in mind that in realistic experiments, 
like atom-optical ones the Hamiltonian of QKR, Eq. 
may be an oversimplification. Therefore, some restric- 
tion on the validity of the present result will be imposed. 
For example, the cold atoms experiment involves colli- 
sions between the atoms, which lead to another dephas- 
ing mechanism. Let us estimate the corresponding col- 
lision dephasing time, r^. The two-particle scattering 
mean free path is known to be Is ~ \/{na?'), where n is 
the atomic concentration and a is the s-wave scattering 
length. The corresponding dimensionless scattering time 
is Tg = Is/ {Tv), where w is a typical atomic velocity that 
may be estimated as f ~ ft/cL |Z|/Sm « HkL \/DciTs/km . 
This leads to the self-consistent estimate of the dimen- 
sionless dephasing time: 



{lskLlKf'\ 



(125) 



Once again to observe the classical-to-quantum 
crossover the inequality 41 e ^ Ts should be valid. 



VIII. CONCLUSIONS 

In this paper, we developed an analytical theory to 
incorporate systematically the Ehrenfest time, Ie, into 
the weak dynamical localization. We map the loop ex- 
pansion, central to the theory of weak dynamical local- 
ization onto the interference of paths in the configura- 
tion space with a certain loop geometry. To propagate 
along such paths a time longer than 2mtE {'m — 2,3, ■ ■ ■) 
(where m is determined by a specific loop geometry) is 
needed. We establish thus that the onset of the dynam- 
ical localization is delayed by the multiples of 2tE- In 
particular, for QKR, the delay is AIe, while for systems 
with broken time-reversal symmetry it is GIe- At shorter 
times, quantum corrections to the linear dispersion (clas- 
sical diffusion) do exist. However, they results only in 
the renormalization of the frequency-independent diffu- 
sion coefficient, and thus are Koi responsible for the onset 
of the dynamical localization. 
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Our quantitative predictions are based on the loop 
expansion and are essentially perturbative. They are 
closely analogous to the weak Anderson localization in 
ballistic systems. They may be considered as a further 
support for the long-standing conjecture of Ref. 13 about 
the similarity between the dynamical localization and 
Anderson localization for irrational k/Ai: The per- 
turbative corrections are responsible only for the early 
evolution of the dispersion function, i.e., the suppression 
of the classical diffusion. At longer time t > t^, the 
strong dynamical localization is expected to develop, as 
supported by arguments based on supersymmetric diffu- 
sive a modeli?. However, it still remains a challenge to 
prove this conjecture in the presence of finite Ie- 

Although our analytical treatment is developed for the 
QKR, we believe that it may be extended to more general 
quantum driven systems. The essential requirements im- 
posed on the dynamics of an underlying classical system 
are: (i) area-preserving; and (ii) the existence of classical 
stochastic diffusion (subject to certain symmetry). Pro- 
vided (i) and (ii) are satisfied, we expect the functional 
form of the weak localization corrections to be the same 
as derived above. The specific of chaotic motion enters 
through modifying the classical quantities such as diffu- 
sion coefficient, Lyapunov exponent and its fluctuations 
- such problem is purely classical and, in general, may 
be resolved (say numerically). 

Technically, the most important part of this work is the 
derivation of the one-loop vertex (Hikami box) without in- 
troduction of a regularization. This allows to analyze ac- 
curately the minimal quantum wave packet. As a result, 
the Ehrenfest time is quantitatively defined (with the log- 
arithmic accuracy) as the time needed to expand an ini- 
tial minimal wave packet up to a macroscopic sizeSLS^. 
Alternatively, the minimal wave packet may be analyzed 
within the Moyal formalism.'^^'^^. A study along this line 
has been reported recently in the context of ballistic su- 
persymmetric cr-model'^*' . This supports the conjecture of 
Ref. I23I that the Ehrenfest time should not depend on the 
regularization, since the later is only intended to mimic 
the effect of quantum diffraction. Indeed, in accord with 
Ref. I23, our current results may be fully reproduced by 
introducing a proper regularization'^'' to the supersym- 
metric a model developed for QKR^^. 

The quantitative predictions made for the in- 
dependent classical-to-quantum crossover in QKR may 
be suitable for experimental verifications in various con- 
texts. In particular, already existing experiments on the 
dynamical localization in the energy growth of ultra- 
cold atomic gases have greatly contributed to under- 
standing of this crossover. Our quantitative predictions 
are expected to be accurate in the asymptotic regime 
k < 1 < K . For their quantitative verification, it is thus 
highly desirable to decrease S down to 0.1 — 1. We also 
pointed out that some periodically driven mesoscopic su- 
perconducting structures may be suitable for realizations 
of the QKR. The dynamical localization and classical- 
to-quantum crossover in these systems are observable by 



monitoring charge fiuctuations of the superconducting is- 
land. It is important to mention that, all realistic exper- 
iments introduce noise and thus a finite dephasing time 
T^p . We have shown here that for an observation of the 
tfi-dependent crossover, the condition AtE ^ r^p must be 
satisfied. 
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APPENDIX A: FINITE TIME CORRELATION 
EFFECTS ON THE SELF-ENERGY 



In this Appendix we study the higher order time corre- 
lation effects, starting from the exact quantum density- 
density correlator, Eq. (|12|l . In particular, we clarify that 
in the semiclassical limit, the higher order corrections, 
namely Eq. (|23|l . (see Ref. Ill|) to the diffusion constant, 
i.e., are found. 

For an arbitrary Tc > 1 Eq. |{T5|) is replaced by: 

{I, n = Si,i, + ^ s (z, n Vo (/", n . (ai) 



Here the self-energy, is given by 

= ^e''^"(;|[/"|r')(/|;7"|/")* (A2) 

k=l 



Note that lo — I'q = I, In — = I", and Ik ^ 1'^. for 
< k < n. Here in order to simplify the notation we de- 
note the matrix elements {lk+\U\l(k+i)+) as J/;^^,;,^^!,^ , 
and similarly for their complex conjugates. These matrix 
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elements may be written explicitly as 



and 



Ui 



iK 



de 



(fc+i)+ 

2tt 



exp 



''(fc+i)+ 
2S 



cos 



Hk+1) + 



{I 



k+ - Hfe+l) + j tl{k+l) + 



ur 



de 



(fe+i)- 



2tt 



exp 



**(/c+l)- 



2h 



(A3) 



~ cos e^t+i)- + ^ {ik^ 



{k+l}~j t>(k+l)- 



To proceed further, we introduce the following quantities: 



rrik = {lk+ ~ Ik- )/Ji, qk = iOk+ - Ok- ) /k ■ (A4) 
Then with the substitution of Eq. I|A3|I we rewrite 



as 



dO{k+i)+ d0(k+i)- 



2tt 



2tt 



exp < inik+i 



(A5) 

^(fe+i)+ + ^(fe+i)- 



2iK . kqk+i . 0(k+i)+ + 

sm sm ^ 

k 2 2 



'(fc+i)- 



lk+ + h 



(fc+l)+ + ^(fe+1) 



+i {mk - mk+i) 



^{k+i)+ + (^{k+i)- 



Furthermore, we insert the Fourier transform: 



exp 



2iK . ^ . kq 
— — sm sm — 
k 2 



n sgn q \ ^ ^^^^^ I c- 



(A6) 



into it with sgn denoting the sign of q. Then Eq. (|A5|I 
is substituted into Eq. (|A2|) . With the sum with respect 
to {lk+ + Ik-) /2 ( A: = 0, 1, 2, • • • ) and the integral with 
respect to {9k+ + 6k-) /2 performed, eventually Eq. (|A2|) 
is reduced to 



2K 

Y.{l,l") ^e"^Jo I — sin%o 



(A7) 



XJn 



T.^- EE En 

r=2 ni ,712,- ■■ „nr "tl ■,m2,--- ,mr-l <30i9l r'- ,gr-l k=l 

kqk-1 



2K 

k sgn Qk-i { 



sm ■ 



^^qk~qk-l,-mk ^nik-rrik-i.-nk sgnqk-i 



Shortly later we will see that the following relations, im- 
plied by the two Kroneck's symbols, 



Qk = qk-i - mk 



(A8) 



nik = rrik-i - Uk sgnqk-i 



(A9) 



are essential for the derivation of higher order corrections 
to the diffusion constanliii. It is in order to emphasize 
that they are exact even at the quantum mechanical level, 
although originally found in the classical context. 

Note that above mg = m„ = 0. So far the discus- 
sions above are formally accurate. We assume now that 
the correlator, 'Do{l, I') does not depend on the center-of- 
mass, namely T>o{l, I') = 'Do{l — V). This is supplemented 
by performing Fourier transformation for 'Do{l — I'): 



-Do {I - n 



^PoMe^'^C-O (AlO) 



and imposing the boundary condition as (resulting from 
finite Tc 



(All) 



with if denoting the Fourier component. Then Fourier 
transforming S(Z — I') leads to 



S(^) = e-Jo('^sin^ 



(A12) 



+E^^"^E EE En 

r=2 "15^0 n2,--- ,nr mi ,m2,--- ,mr-i qi,--- ,q, — 2 fc=l 

(2K kip\ (2K kip 
xJ„J— sin- J„ (— sin- 



X J, 



Uk sgnqk-i 



2K 



. kqk-i 
sm 



Note that above we suppressed the two Kroneck's sym- 
bols to simplify the expression. One should keep in mind 
that the sum over m's, g's and n's are restricted by the 
two "motion" equations, i.e., Eq. (|A8|) and ljA9|l . 

Let us make the semiclassical approximation, i.e., 
\kqk\ <C 1 in Eq. HA12p and focus on the limit Kip <C 1. 
The first term in the self-energy leads to the diffusion 
constant as K'^/4:, as discussed in Sec. IIII Bl The sec- 
ond term gives higher order oscillatory corrections. In 
fact, up to (Kip)'^, the sum is contributed by a particu- 
lar series of {qk, mk) {k — 0, 1, 2, • • • ) (so-called Fourier 
paths)" as {(p, 0) ±(0, 1) • • -±(1, -1) {ip, 0). Here 
• • • is shorthand of product of Bessel functions, which is 
an expansion in powers of K~^^^. Since for Kip ^ 1, 
Jn{Kip) ~ (Kip)"', therefore, ni_r = ±1 is the only con- 
tribution to the order (Kip)'^. This then leads to S((p) 



^{^)^e''^MKip) 



(A13) 



+ EE E En 

r=2 rii ,nr = ±l n2,-- - ,"r-l "il,™2,-- - ,™r-l 91 ,9i--2 k=l 

X Jr,i{K(p)JnAK(p)Jnk sgnqk-i m\qk-i\) 
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In the limit ujTc <C 1 , Kip <C 1, one recovers the difFuson 
Eq. 122|l with the diffusion constant given by Eq. (|23|) up 
to 0(1) (for K > 1)". 

Proceeding along this line, we may reproduce Shep- 
elansky's results for quantum diffusion constant of early 
evolution^i*^. The basic observation is that if number 
of kicks is less than 4, then qk = 1. Based on Eq. ljA12p . 
this then implies that the results of Eq. (|23|l still holds 
except that K is replaced by 2^ sin(^/2)/^. 



APPENDIX B: TWO RELATIONS RESULTING 
FROM TRS 

In this Appendix we show two exact relations reflecting 
TRS. In the first, for any n > , 

Po {I, 0- 1', 0'; n) = Vo {I', -9'- 1, -6; n) . (Bl) 

Proof. Use the mathematical deduction. For n — 1 , 

Vo{i,e-i\e'-i) (B2) 

= PS {I -I') 5 [9-9' -I') 

= S{1- Ksm{9 5{9-l-9' - I') . 

On the other hand, 

Vo{l',-9';l,-e;l) (B3) 

= ~P 5{l' ^l) 5{-9' + 9-1) 

= 5{l' + Ksm{9' + I') -I) 5 [-9' - I' + 9 - I) 

= 5{l' + Ksm{9 -l)-l) 6{9 -1-9' -I') . 



Comparing the last two lines of Eqs. ljB2|l and ljB3|) . we 
immediately see that Eq. (jBlll holds for rt = 1 . Next 
we assume that Eq. (jBlll holds for arbitrary n = k > 1 . 
Then for ri = fc + 1 , we obtain: 

Va{l',-9';l,-9;k+l) (B4) 

= PVoil',-9';l,-9;k) 

= PVo{l,9;l',9';k) 
= Voil,9;l',9';k + l) . 

Thus, Eq. (inU also holds for n = fc + 1. Q. E. D. 
Now we turn to show the other relation. That is, if 



then 



Proof. 



fil,9;l',9') = f{l',-9';l,-9) , (B5) 

P f{l,9:l',9')^ f{l,9;l',9')PT. (B6) 

f(l,9; I', 9') Pt 

fil,9; I' + Ksm{9' + l'),9' + 1') 

f {I' + K shi{9' + I'), - {9' + I'): 1,-9) 

Pfil',-9'; 1,-9) 

P f{l,9-l',9') , (B7) 



where we used the definition, Eq. H30|l . in the second line, 
Eq. (|B5() in the third and fifth lines, and the definition, 
Eq. lini), in the fourth line. Q. E. D. 



APPENDIX C: DERIVATION OF 
RENORMALIZATION FACTOR OF DIFFUSIVE 
COOPERON 

In this Appendix, we derive the renormalization fac- 
tor, Eq. H37|) . describing modification of the diffusive 
Cooperon due to the propagation through the Lyapunov 
region. We first analyze the asymptotic instability of a 
generic chaotic trajectory. 



1. Asymptotic instability 

By varying the equations of motion of the classical 
kicked rotor: 



d9 
'dt 
dl_ 
dt 



Ks\n9^5{t-n) 



and defining the variables: 
obtain: 



\ii\59\, a 



dz 
'dt 



da 



— -l-a^ = Kcos9y^S{t- n), 



(CI) 
61/69, we 

(C2) 



which describes the evolution of separation of two nearby 
trajectories along a reference trajectory, initiated from 
{Iq , 9o). From Eq. ljC2p we see that a is a fast changing 
variable. That is, the dynamics of a introduces some clas- 
sical time scale, beyond which the slow changing variable 
- z is independent of initial a. Let us further introduce 
the variables: ari~denoting a right after the n-th kicking, 
and z„-denoting z at the n-th kicking. Equivalently, we 
rewrite Eq. HC2|) as 



Zn+l 



Zn = ln(l + an) , 
a:\n+i 



a 



1 



K cos 9, 



n+l ■ 



(C3) 



For n » 1, (z„) = An, ((z„ — Xn)^) = 2X2 n. Thus, the 
Lyapunov exponent A and its fluctuation A2, character- 
izing the long time instability are defined as 



A — lim A(n), A2 = lim A2(n), 

n — >oo n — >oo 

n — 1 

^(") = - E In (! + «!„') , (C4) 



A2(n) 



n'=0 

E 

.n'=0 



In (1 + a|„') — ?iA 
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respectively. Note that at finite times, they are 
trajectory-dependent. In the Hmit n — > oo, the time av- 
erage is expected to be equivalent to the average over a, 
as well as the initial conditions located in the stochastic 
region. 

The estimation of A and A2 in the limit K ^ \ may 
be made as follows based on the simple analysis above: 



= HK/2) ■ (C5) 
A2 = (In^ |i4:cos6i|) - A^ = C(3) -ln^2 w 0.82, 



A = (ln|ifcos( 

,2 



where the angular brackets imply uniform averaging over 
the angle. 



2. Renormalization factor 

According to the definition of 71^ , Eq. (|34|l it is the kick 
number when z « with the logarithmic accuracy. That 
is. 



ric-l 

= zo + X! 1^(1 + 

i=o 



(C6) 



equivalently. 



ric — 1 

Z0+ 51 [ln(l + a|j)- A(nc)] 



Here zq is some function of the initial deviations 
((5^0, (5^0 ) • The exact form of zq is unessential. For 
an estimate, we notice that in the limit K ^ 1 
the evolution of 59n may be approximated as {50i^ + 
61q/ K cosOq) Hfe^o KcosOj^. Therefore, 

Zq = In I (5^0 + (^'o/-^ cos ^0 1 

~ In ^5el + 5ll/K^. (C8) 



Substituting Eq. (|C7|I into exp(2ici;nc), we obtain: 
2iuj zq 



X exp 



A (71c) 

2iw 



(C9) 



n-c — 1 



A(nc) 



[ln(l + a|,)-AK)] 



3=0 





2iw Zq 




exp 


A(nc) 


exp 







A^K) 



where in the second line we use the fact that 
ujX2{nc)nc/ X^iric) ^ 1 . In the limit ^ 1 , A(nc) 
A, A2(nc) A2 . Moreover, from Eq. (jC7|) we see that 
nc^t'^ = -zo/A . Thus, Eq. ifUO)! is reduced to Eq. (|S7jl . 



(a) 





(b) 




(c) 




FIG. 9; Diagrams that lead to the interaction vertex with 
the particle conservation law respected: (a) (JOi, (b) 5'D2, and 
(c) SV-i. 



APPENDIX D: THE EXACT INTERACTION 
VERTEX 

In this Appendix we discuss the exact interaction ver- 
tex appearing in the one-loop calculation. The general 
diagram, sketched in Fig.|21 may be categorized into three 
classes, Fig.^a-c. In order to make the formula compact, 
let us write them down in the operator representation as: 



51)2 = VvCT), 
ST>3 = VGVv- 
We then insert the identity : 

e*" -P = 1 - 



'V -I 



into Eq. HD1() and rewrite SVi as 
S'Di='DC'D + 5V[ 

with 



(Dl) 

(D2) 
(D3) 

(D4) 
(D5) 



6t>[ = 

5Vi = 



1 V 



CV + VC 



1 V 



V -1] V C 



'V-1] V 



(D6) 



For the first term in Eq. (jPSp , it was proved by Altlandi^ 
that it vanishes in the quantum limit. In the next section, 
we show that this term is indeed pure classical, and docs 
not contribute to the quantum interference correction. 
By the same token we show that the loop expansion does 
not violate the particle conservation law. That is the 
form of the interaction vertex implies that the quantum 
corrections may be reduced to the renormalization of the 
diffusion coefficient. 
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1. Semiclassical analysis on VC'D 



A cancellation mechanism 



Under appropriate approximation, it was proved that 
2? C I? is included into classical dijfusive propagatOB^S. 
Thus it does not violate diffusion equation. Since we 
are interested in the dynamics involving Ehrenfest time, 
the diffusive propagator can not serve as a starting point 
of the formalism developed here. A natural question is 
whether this conclusion may still be applicable. Although 
it remains a challenge to prove at the accurate level, here 
we present a physical interpretation. We conclude that 
this term is pure classical, included in the exact classical 
propagator-solution of the FPR equation. 

For quantitative discussions below, let us denote 
the angular momenta and angles, appearing in the 
retarded/advanced propagator lines of Cooperon as 
h±,l2±r ■■ , ln± and 6'i± , 6'2± , • • • , 6'(„-i)± following 
the forward time direction of the retarded line (see 
Fig. |31 for notations). Counting from the most left 
hand side of the Cooperon and following the backward 
time direction of the retarded line, we denote those as 
^o± , l-i± , l~2± , • • • , and 9o± , 9-i± , 9_2± , • • ■ ■ Semi- 
classically, only the constraint below 



\0 



1+ + C'(„_i)_ 



1^: 



2+ + C'(„-2)- 



l^(n-l 



(D7) 
is imposed. 

According to the exact Eq. ljA8|l . at the boundary of 
Cooperon and (left) Diffuson, 

{01+ - - {00+ - 00-) - {li+ - 01-) . (D8) 

On the other hand, Eq. ljD7|l does not impose any con- 
straints on 01+ — 01-. If it is order the usual Wigner 
transform may be performed (see Eq. (|14|l '). That is to 
say, the left-most kicks of the Cooperon may be incorpo- 
rated into the left Diffuson. 

This discussions remain applicable for all the suc- 
cessive pairs of kicks of the Cooperon, until we meet 
some fc > 1, where 0k+ ~ 0k- ^ 1-a typical feature 
of the Cooperon. This is equivalent to the validity of 
Eq. (|D8p with the right hand side of order 1. In fact, 
(^0+ — So ) and {li+ — 0i-) develop from the left (classi- 
cal) Diffuson-foUowing Eqs. HA8|I and HA9|I . In the Diffu- 
son side, {0k+ - 0k-) {k = -1, -2, • • • ), {lk+ - h-) {k = 
0,— 1,— 2,---) are of order k. Moreover, any evolution 
can last at most finite classical correlation time Tc- Thus, 
(6*1+ — 01-) must be of order k. Thus, Eq. IjDSp can not 
be satisfied ff {0i+ - 0i-) - 0(1). 

Thus we may conclude that I? CP is pure classical. 
It characterizes the classical probability of trajectories 
with a peculiar feature. Since \0i- + 0n+\ ^ k <s:_ 1, 
01+ ^ 01- ^ —0n+, the trajectory switches its angle to 
the direction (almost) opposite to the initial. 



Having the general expression, Eq. H52|l at hand, we 
turn to show that SI)'^ , S'D2, SV^ cancel each other at 
the semiclassical level, namely: 

5V[ + 5V2 + (5X>3 = . (D9) 

Indeed, employing identity, Eq. ljB6|) . we obtain: 

5V[ {l,0;l\0') 



vlVoi l,0;l" + 



2 ' 



5V2 {l,0;l\0') 



-Do 



2 ' 



-0" 



(DIO) 



(Dll) 



|27rS(5 


I - 




S 











r' + '-^\-i 



^Vo[l"-^-^,-0"+^-^;l',0' 
u 1 2 ' 2 



and 



SV3 {l,0;l',0') = V 



5I2 



Vo[lAl"+'-^,0"+^-^ 



x2TrkS \l" - -^-n 6 - 



— - 6*' - - — 
2 2 
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Since Vq is the solution of FPR equation, i.e., Eq. l(TH|l . 
these three terms cancel each other. Therefore we find 
that only SV^ leads to the non-vanishing one-loop cor- 
rection to the classical density-density correlator. 



APPENDIX E: ALEINER-LARKIN 
REGULARIZATION IN KICKED ROTOR 

It is easy to see that in the classical limit: S — > 0, 
the functional form of SV^ and (52?2,3, i-e., Eqs. (jD10|l . 
(|D11|I and (|D12|I are identical to what have been found 
for classical Lorentz gas (leaving the feature of standard 
map aside). These terms cancel each other, leading to 
the absence of the weak localization correction. Remark- 
ably, this cancellation still holds, even if the initial min- 
imal wave packet is taken into account. On the other 
hand, the weak localization correction does exist, given 
by the exact vertex SV^, which is absent in the expansion 
of ballistic supersymmetric a model (without a regular- 
izer). In this way, one may wonder that an appropriate 
regularization may lead to a physical description of weak 
localization in semiclassical chaotic systems. This, in- 
deed, was developed by Aleiner and Larkin for ballistic 
electronic problematic. However, an important issue 
remains open whether and to what extent the physical 
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results depend on the regularizer, rather than on the in- 
trinsic quantum nature of the problem. The exact in- 
teraction vertex may serve as a testing ground of this 
regularization. 

In fact, the regularization introduced in Refs. I27i58l is 
also applicable for QKR. To this end we try to mimic 
the "Born impurity" by modifying the (quantum) free 
rotation operator Vv in the foUowing way: 



with the regularized Diffuson V and the Cooperon C sat- 
isfying the following equation: 



Vv Vv exp 



i6S 



(El) 



Here 5S is some stationary random perturbation that 
commutes with Vv i i-e., [5S, Vv] = 0. Moreover, we 
assume that it is short-ranged correlated in the angular 
momentum space: 



1 - 



1^ 
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{did Si dv 5Sv) ^-5i,v 
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In the classical limit, this additional stationary random 
perturbations leads to the modification of the standard 
mapping in the following way 



In+l = In + KsinOn, 



9;„+i 5Si„ 



(E3) 



Then following the same procedure as in Ref. 27, we find 
the one-loop correction to the density-density correlator 
to be 



5v {i,e;i',e') 

2 r dli dO 

T, 



-C (Zi, 6*1; Zi, - 



(E4) 



d 
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Eqs. ljE4p and (|E5p are fully analogous to those found 
for ballistic electronic systems. In Eq. (|E5(I . the regu- 
larizer of FRP equation mimics the spread of minimal 
wave packet arising from intrinsic quantum diffractions. 
Originally, it was expected^^ that the only physical effect 
arising from this regularizer is to determine the Ehren- 
fest time since it smears the sharp classical propaga- 
tor. Indeed, one may further follow the procedure of 
Ref. I23 to calculate the weak-localization correction to 
the diffusion constant. As a result, the functional form 
of Eq. (|82|l is reproduced with Ie acquiring explicit Tq- 
dependenceSl. At (Ar,)^^ ^ ^/K, the result thereby 
obtained are identical to Eq. This reflects a basic 

belief previously anticipated^LSS^i. That is, a physical 
regularizer strength must match the minimal quantum 
wave packet. 



^ D. M. Basko, M. A. Skvortsov, and V. E. Kravtsov, 

Phys. Rev. Lett. 90, 096801 (2003); V. E. Kravtsov, 

cond-mat/0312316 
2 F. L. Moore et al, Phys. Rev. Lett. 75, 4598 (1995); C. F. 

Bharucha et al., Phys. Rev. E 60,3881 (1999). 
^ M. G. Raizen, Adv. At., MoL, Opt. Phys. 41, 43 (1999). 
* C. Zhang, J. Liu, M. Raizen, and Q. Niu, Phys. Rev. Lett. 

92, 054101 (2004). 
^ O. A. Starykh, P. R. J. Jacquod, E. E. Narimanov, and A. 

D. Stone, Phys. Rev. E 62, 2078 (2000). 
^ G. Casati, B. V. Chirikov, J. Ford, and F. M. Izrailev, 

in Stochastic Behavior of Classical and Quantum Hamilto- 

man Systems, Lecture Notes in Physics 93, edited by G. 

Casati and J. Ford, 334 (Springer, New York 1979). 

B. V. Chirikov, F. M. Izrailev, and D. L. Shepelyansky, 

Sov. Sci. Rev. Sec. C 2, 209 (1981). 
® H. Ammann et al., Phys. Rev. Lett. 80, 4111 (1998). 
® B. V. Chirikov, Phys. Rep. 52, 263 (1979). 
° A. L. Lichtenberg and M. A. Lieberman, Regular and 

Chaotic Dynamics (Springer-Verlag, Berlin 1991). 



" A. B. Rechester and R. B. White, Phys. Rev. Lett. 44, 
1586 (1980); A. B. Rechester, M. N. Rosenbluth, and R. 
B. White, Phys. Rev. A 23, 2664 (1981). 
M. Khodas and S. Fishman, Phys. Rev. Lett. 84, 2837 
(2000); Erratum, tbid. 84, 5918 (2000); M. Khodas, S. 
Fishman, and O. Agam, Phys. Rev. E 62, 4769 (2000). 
S. Fishman, D. R. Grempel, and R. E. Prange, Phys. Rev. 
Lett. 49, 509 (1982). 

S. Fishman, D. R. Grempel, and R. E. Prange, Phys. Rev. 
A 29, 1639 (1984). 

R. Bliimel and U. Smilansky, Phys. Rev. Lett. 69, 217 
(1992). 

P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 57, 
287 (1985). 

L. P. Gorkov, A. I. Larkin, and D. E. Khmelnitskii, Pis'ma 
Zh. Eksp. Teor. Fiz. 30, 248 (1979) [JETP Lett. 30 , 248 
(1979)]. 

A. Altland, Phys. Rev. Lett. 71, 69 (1993). 

A. Altland and M. R. Zirnbauer, Phys. Rev. Lett. 77, 4536 

(1996). 



27 



K. Efetov, Supersymmetry in disorder and chaos (Cam- 
bridge University Press, UK 1997). 

A. I. Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 
55, 2262 (1968) [Sov. Phys. JETP 28, 1200 (1969)]. 

22 G. P. Berman and G. M. Zaslavsky, Physica A 91, 450 
(1978); G. M. Zaslavsky, Phys. Rep. 80, 157 (1981). 

23 For a review, see F. M. Izrailev, Phys. Rep. 196, 299 
(1990). 

2* Ph. Jacquod and E. V. Sukhorukov, Phys. Rev. Lett. 

92, 116801 (2004); J. Tworzydlo, A. Tajic, and C. 

W. J. Beenakker, Phys. Rev. B 69, 165318 (2004), 

cond-mat/0405122 
2^ M. C. Goorden, Ph. Jacquod, and C. W. J. Beenakker, 

Phys. Rev. B 68, 220501 (R) (2003). 
2*^ J. Tworzydlo, A. Tajic, H. Schomerus, and C. W. J. 

Beenakker, Phys. Rev. B 68, 115313 (2003); P. G. Sil- 

vestrov, M. C. Goorden, and C. W. J. Beenakker, Phys. 

Rev. Lett. 90, 116801 (2003). 
2^" L L Aleiner and A. L Larkin, Phys. Rev. B 54, 14423 

(1996). 

28 M. G. Vavilov and A. I. Larkin, Phys. Rev. B 67, 115335 
(2003). 

2® C. Tian, A. Kamenev, and A. Larkin, Phys. Rev. Lett. 93, 

124101 (2004). 
3° J. Miiler and A. Altland, nlin.CD/0411055 
3^ D. L. Shepelyansky, Physica D 28, 103 (1987). 

32 D. A. Steck, V. Milner, W. H. Oskay, and M. G. Raizen, 
Phys. Rev. E 62, 3461 (2000). 

33 P. H. Jones et. al, physics/0405046 P. H. Jones et. al, 
quant-ph/0309149 

3^ S. Montangero, A. Romito, G. Benenti, and R. Fazio, 

cond-mat/0407274 
3^^ A. Isacsson, L. Y. Gorehk, R. L Shekhter, Y. M. Galperin, 

and M. Jonson, Phys. Rev. Lett. 89, 277002 (2002); L. Y. 

Gorehk, A. Isacsson, Y. M. Galperin, R. L Shekhter, and 

M. Jonson, Nature 411,454 (2001). 
38 S. Hikami, Phys. Rev. B 24, 2671 (1981). 
^"^ M. R. Zirnbauer, in Supersymmetry and Trace Formulae, 

Chaos and Disorder, edited by I. V. Lerner, J. P. Keating 

and D. E. Khmelnitskii (Kluwer Academic/Plenum, New 

York 1999). 

38 See, e.g., G. M. Zaslavsky, M. Edelman, and B. A. Niyazov, 

Chaos 7, 159 (1997). 
3^ F. M. Izrailev, Phys. Rev. Lett. 56, 541 (1986). 
*° M. Thaha, R. Bliimel, and U. Smilansky, Phys. Rev. E 48, 

1764 (1993). 



T. Tworzydlo, A. Tajic, and C. W. J. Beenakker, 
cond-mat/0405122 

"'2 R. A. Smith, I. V. Lerner, and B. L. Altshuler, Phys. Rev. 
B 58, 10343 (1998). 

"3 R. S. Whitney, I. V. Lerner, and R. A. Smith, Waves Ran- 
dom Media 9, 179 (1999). 

C. Tian and A. I. Larkin, Phys. Rev. B 70, 035305 (2004). 
For a review see e.g. S. Chu, Science 253, 861 (1991); W. D. 
Phillips, Rev. Mod. Phys. 70, 721 (1998); C. E. Wieman, 

D. E. Pritchard, and D. J. Wineland, Rev. Mod. Phys. 71, 
S253 (1999). 

^8 R. Graham, M. Schlautmann, and P. ZoUer, Phys. Rev. A 
45, R19 (1992). 

F. L. Moore et al, Phys. Rev. Lett. 73, 2974 (1994); J. C. 
Robinson et al, Phys. Rev. Lett. 74, 3963 (1995). 
"8 In the conventional QKR is periodic, i.e. *(6l + 27r) = '^{9), 
leading to the discrete angular momentum: I = nk with n 
integer. Periodically kicked cold atomic gases differs from 
it, with the boundary condition replaced by: ^'(6 -I- 2tt) = 



*(6l)e 



where lo G [0, 1] is a fractional part of an atom's 



momentum in units 2hkL- It may be eliminated by the 
gauge transformation that introduces an Aharonov-Bohm 
flux into the kinetic energy {l—lo)^/2. Such flux does not af- 
fect dynamical localization, though quantum resonances2£ 
is eliminated by averaging over Iq. 

C. Tian, A. Kamenev, and A. I. Larkin, Bull. Am. Phys. 
Soc. 49, 482 (2004). 
^° E. Chow, P. Delsing, and D. B. Haviland, Phys. Rev. Lett. 
81, 204 (1998). 

^'^ D. B. Haviland and P. Delsing, Phys. Rev. B 54, R6857 

(1996) . 

^2 M. Watanabe and D. B. Haviland, in Studies of High Tem- 
perature Superconductors (Advances in Research and Ap- 
plications), Vol. 43, edited by A. Narlikar (Nova Science 
Publishers, New York, 2002). 

^^3 D. Cohen, Phys. Rev. A 44, 2292 (1991). 

"^'^ V. Milner et. al, Phys. Rev. E 61, 7223 (2000). 

'^^ B. G. Klappauf et. al, Phys. Rev. Lett. 81, 1203 (1998). 

''^ A. J. Daley and A. S. Parkins, Phys. Rev. E 66, 056210 
(2002); G. Duffy et. al, cond- mat/0401346 

^"^ K. B. Efetov, G. Schwiete, and K. Takahashi, Phys. Rev. 
Lett. 92, 026807 (2004). 

'^^ I. L. Aleiner and A. I. Larkin, Phys. Rev. E 55, R1243 

(1997) . 



